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ABSTRACT 



The Far Ultraviolet Spectroscopic Explorer [FUSE) has allowed precise deter- 
minations of the column densities of molecular hydrogen (H2) in Galactic lines of 
sight with a wide range of pathlengths and extinction properties. However, sur- 
vey studies of lines of sight with greater extinction have been mostly restricted 
to the low-J states (lower total angular momentum) in which most molecular 
hydrogen is observed. This paper presents a survey of column densities for the 
molecular hydrogen in states of greater rotational excitation (J > 2) in Galactic 
lines of sight with logA^(H2) > 20. This study is comprehensive through the 
highest excited state detectable in each line of sight. J = 5 is observed in every 
line of sight, and we detect J = 7 in four lines of sight, J = 8 in one line of sight, 
and vibrationally excited H2 in two lines of sight. We compared the apparent 
6- values and velocity offsets of the higher- J states relative to the dominant low-J 
states and we found no evidence of any trends that might provide insight into 
the formation of higher- J H2, although these results are the most affected by the 
limits of the FUSE resolution. We also derive excitation temperatures based on 
the column densities of the different states. We confirm that at least two distinct 
temperatures are necessary to adequately describe these hues of sight, and that 
more temperatures are probably necessary. Total H2 column density is known to 
be correlated with other molecules; we explore if correlations vary as a function 
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of J for several molecules, most importantly CH and CH+. Finally, we briefly 
discuss interpretations of selected lines of sight by comparing them to models 
computed using the Meudon PDR code. 



INTRODUCTION AND BACKGROUND 



H2 is the most common molecule in the Universe, as well as a signiflcant diagnostic of 
dense and dusty conditions in interstellar clouds, critical to the physics and chemistry of star- 
forming regions. The many rotational-vibrational-electronic transitions of H2 are found in 
the far-ultraviolet (FUV; A < 1200 A) portion of the spectrum. Thus, it was not until the in- 
troduction of space-based instruments with sufficient sensitivity and spectral resolving power 
in the FUV that detailed observations of H2 were undertaken. While instruments such as 
Cope rnicus were capable of measuring H2 (e.g. ISpitzer et al.lll973l Il974j : ISpitzer fc Cochran 



19731 ). observations were generally restricted to lines of sight with Ay ^ 1 mag. The most 
recent mantle-bearer of FUV — and therefore H2 — observations was the Far Ultraviolet Spec- 
troscopic Explorer {FUSE), which operated from 1999 through 2007. FUSE was able to push 
well past the extinction limits of previous FUV instruments, observing many Galactic lines 
of sight with Ay ~ 2 — 3 and even a few targets with Ay ~ 5. These observations have 
established a somewhat surprising ubiquity of H2, as it is found to be not just in dense, 
star-forming regions but is also often present in signiflcant am ounts throughout the diffuse 



ISM, both in the Galactic disk and at higher Galactic latitudes ( lWakkerll2006l : iGillmon et al. 
20061 . and references therein). 



Rotational, vibrational, and electronic splitting occurs in the energy levels of H2. The 
transitions observed in the FUV are dominated by H2 that is in the ground electronic state, 
the ground vibrational state, and states of low rotational energy. However, a detectable 
amount of rotationally excited, or "high-J"0 H2 typically exists in the lines of sight that 
have been observed with FUSE. 

The vast majority of the Copernicus and FUSE studies mentioned above measured 
either H2 in all detected J levels for lower column density lines of sight in the diffuse ISM 
(iV(H2) < 10^°) or only J = - 1 H2 in denser lines of sight (A^(H2) > 10^°). The recent 



survey of ISheffer et al.l ( l2008l ) explored H2 in 58 lines of sight using FUSE data, but the 
emphasis of that survey was on the chemical relationsh ips between Ho and other molecules, 
including ^^CO, CH, CH+, and CN. In their paper, ISheffer et al.l ( l2008l ) only indirectly 



More technically, J is the quantum number of the total angular momentum — rotational and electronic. 
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reported H2 column densities as a function of J level by reporting excitation temperatures 
up to To4, from which the individual colu mn densities can be calculated (but only up to 
J = 4). Furthermore, ISheffer et al.l (|2008[ ) did not address issues specific to the excitation 
of the various J levels. 

In contrast, the aim of this paper is to carry out a survey inclusive of all detectable J 
levels and any detectable vibrationally excited material. (In our survey this turns out to be at 
least J = 5 in every line of sight; see §3] and Table IH) We will do this by exar nining reddened 



lines of sight frorn the ■FL/.S'.E archives, primarily from a serie s of papers (IRachford et al. 



200ll . l2002bl I2OO9I : ISnow et al.ll2000l : ISonnentrucker et al.ll2002[ ). many of which were led by 
authors on this paper (BLR and TPS), that used FUSE to examine lines of sight thought to 
potentially contain translucent c louds. We will hereafter ref er to this collection of papers, 
with emphasis on the two papers ( IRachford et al.ll2002bl . 120091 ) that were surveys of J = — 1 
H2 in these lines of sight, as the FUSE Translucent Cloud Survey, or FTCS. Our intent 
is to supplement the measurements of A^(0) and A^(l)@ in the FTCS survey papers with 
measurements of A^(> 2). We will loosely call H2 in the J > 2 states "rotationally excited" 
or "high- J" at points in this paper, although a better physical distinction occurs at about 
J = 3 (see §4.2p . The target sample is discussed more in ^ 

There are several issues that our study intends to address. First, any study based on 
a data set of reddened lines of sight may help better characterize the nature of translucent 
clouds. Qualitatively, translucent clouds are a transition phase between diffuse and dense 
clouds in the ISM. In the diffuse ISM, the interstellar material is dominated by neutral atoms 
and ions, hydrogen is nearly all atomic, dust is present but gas dominates the material, and 
ices do not form on grains. In contrast, in the dense ISM, clouds are dominated by molecules 
(including molecular hydrogen), the depletion of atoms directly onto dust grains is significant, 
and ices readily form on grains. In translucent clouds, a transition between these properties is 
expected to be observed — atoms and molecules will have comparable abundances and chemi- 
cal significance, the molecular fraction of hydrogen [/(H2) = 2A^(H2)/(A^(H I)+2A^(Il2))] will 
approach unity, and ice formation on dust grains begins. In addition to the F TCS papers and 
references therein, recent discussion of translucent clouds can be found in Snow &: McCall 



( I2OO6I ): in particular, note Fig. 1 of that paper. 



In this context of translucent lines of sight, we can further address simple first-order 
questions about highe r- J H2 in the ISM. For example, preliminary work on this topic 
(IRachford et al.l l2002al ) suggested that higher- J H2 might exist in a line of sight primar- 



^Our notation from this point forward is that N{J) refers to the column density of H2 in the J state (or 
states if a range of J is specified). The alternate, condensed notation Nj is used in our figures and tables. 
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ily in a warmer shell around the colder cloud core containing the lower- J H2. Another 
simplistic possibility is that there two distinct clouds along a line of sight, one cold and one 
warm. The hypothesis that higher- and lower- J H2 are not well- mixed can be tested in part 
by analyzing the 6-values (the velocity dispersion if a single distribution is assumed) and 
velocity offsets of the different J levels — differences in either or both would provide circum- 
stantial evidence that the low-J and high-J H2 is, at least in part, physically distinct. For 
example, if there is detectable, systematic offset in velocity offsets between high- and low-J 
H2 in a given line of sight, that might indicate a situation similar to or approximated by 
the two-cloud case. A distribution of such offsets or differing 6-values over many lines of 
sight that cannot be explained as statistical variation could indicate the degree to which the 
different J- levels are mixed. Note that a null result, however, does not necessarily indicate 
the converse, that high- and low-J H2 are well-mixed in most or all cases. 

F urthermore, the ques tion of how higher- J H2 forms and is sustained in the ISM remains 



open. iNehme et al.l (120081 ) examined the FUSE sp ectrum of HP 1020 6 5, and attemp ted to 
model the line of sight with the Meudon PDR code (ILe Petit et al.ll2006l ). iNehme et al.l found 
that the Meudon PDR code could reproduce some of the basic observational properties 
of this line of sight, such as the low-J column densities (including the ratio of ortho- to 
parahydrogen) and the column densities and emissivities of certain molecules. However, 
there are other properties that are not explained, such as the column densities of J > 3 H2 
and CH+. 

The scope of this paper is to present measurements of A^(> 2) in 22 selected lines of 
sight from the FTCS. We will also present some interpretations and simple modeling of our 
results, using the Meudon PDR code. 



OBSERVATIONS AND DATA REDUCTION 



The spectra under examination are taken from the FUSE archives. As noted above, 
our primary source for lines of sight under consideration was the FTCS. This comprises a 
database of about 40 lines of sight. From this database, we have selected a sample of 22 lines 
of sight. Our first selection criterion was a signal-to-noise ratio (S/N) large enough to detect 
the weaker lines of H2 in the J = 4 — 5 states, combined with a qualitative assessment of the 
continuum that would allow for the simple measurement of most lines. Our second selection 
criterion was that the line of sight's H2 content had not been previously examined in detail. 
This latter criterion eliminated five lines of sight (of the original 40) from consideration — HD 
73882, HD 96675, HD 102065, HD 108927, and HP 110432. HD 73882 and HD 110432 were 
analyzed as part of the FTCS JSnow et al.l[200ol : iRachford et al.ll200ll . respectively), while 
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the other three hnes of sight were analyzed in iGry et aLl ( 120021 ) . Furthermore, four of these 
hues of sight (all but H P 73882) have had their column densities used as the sub j ects o f 
modeling — HD 102065 in lNehme et al.l ( 120081 ) and the other three in lBrowning et al.l ( 120031 ) . 



Of the lin es of sight that we d o include, a few have b een analyzed elsewhere, no tably 



HD 185418 in 



Sonnentrucker et al.l (120031 ). HD 192639 in ISonnentrucker et al.l (120021) and 



Lacour et all (l2005|), and HD 206267 and HD 207538 also in llacour et al.l tood ). How- 



ever, we have reanalyzed these lines of sight f or various reasons. Our fitting method was 
d ifferent than ISonnentrucker et al.l (120021 . |2003|) and we detect J = 6 in HD 192639 whereas 



Sonnentrucker et al. 



J2OO2I ) do not report a detection. Furthermore, the analysis of the higher 
J levels is not discussed to the same degree in these papers as it is in the pap ers examining 
the ex cluded lines of sight mentioned above. The reason for reanalyzing the iLacour et al. 



( 120051 ) lines of sight is that our methodology is different as we are using additional curve- 
of-growth methods to explore various possibilities for the velocity structure. The lines of 
sight that are included in our study are shown in Table [1] along with some basic line of sight 
parameters. 



2.1. FUSE Bata 



The Ft/^E" satellite recorded data on four different physical detectors (LiFl, LiF2, SiCl, 
and SiC2), named after the two different reflecting materials on the optical elements — lithium 
fluoride and silicon carbide. Each detector has two adjacent segments, covering different 
wavelengths, denoted A and B (e.g. LiFl A and LiFlB). Fo r fu rther informat i on ori FUSE 
instrumental details and performance, see iMoos et al.l ( I2OOOI ) or ISahnow et al.l (120001 ). 



The spectra were processed in the CALFUSE pipeline, using version 2.4.0 or later. All 
of the analysis described below was conducted in different detector segments independently, 
before merging results from various segments. This was done to avoid the possible effects of 
any i nstrumental features because, unlike a study using a limited number of absorption lines 
(e.g. I Jensen &: Snowll2007l ). the lines in this study cover nearly the entire FUSE wavelength 



region. Therefore, a global coadd of different detector segments is sure to include at least 
some detector segment-specific instrumental features. Analyzing the detector segments in- 
dividually, therefore, provides us with the ability to look note possible instrumental effects. 
Some minor unidentified features were noted that influenced whether or not some absorption 
lines were included, but this affected, at most, one or two lines per J level. 
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2.2. Auxiliary Data 

In order to interpret our results, it is useful to have several pieces of auxiliary informa- 
tion, including but not limited to A^(0) and A^(l). These low- J column densities of H2 are 
taken from the FTCS, except for HD 195965, the only line of sight in our sample which is 
not in the FTCS. For this line of sight, A^(0) and A^(l) are taken from a FUSE survey of H2 
in the Galactic disk (J. M. ShuU 2010, in preparation). Additional observational values for 
these lines of sight, such as molecular column densities and values of reddening and extinc- 
tion, are taken from various sources. Reddening and extinction values, including references, 
are given in Table [21 

There are literally hundreds of H2 transitions which have wavelengths corresponding to 
the spectral coverage of FUSE. In this study, we typical measure ~ 80 lines of J > 2 H2 per 
line of sight (though more may be observed; see ^2.31). Wa v elength s and oscillator strengths 



for the H2 absorption lines were taken from I Ab grail et al.l ( ll993al Jbl). These data were used 



to identify observed transitions, derive velocity offsets, and make curve-of-growth fits (see 



2.3. Fitting of H2 Lines 



We follow our previously established methods, most recently outlined in I Jensen Sz Snow 



( 120071 ). Our goal is to measure equivalent widths for use in curve-of-growth determinations 
of column densities. This is detailed below in §2.51 To summarize, we measure the equivalent 
widths of observed absorption lines by making simple fits to the local continuum and struc- 
ture of the lines. The local continuum around each line is fit with a low-order polynomial. 
Weaker absorption lines, where the observed profile is dominated by the instrumental profile, 
are fit with Gaussians. Absorption lines where the observed profile depth approaches unity 
are fit with generalized Voigt profiles, assuming a 15 km s~^ full-width at half-maximum 
(FWHM) Gaussian instrumental profile for FUSE. Errors are derived from standard error 
propagation in the Gaussian case, including continuum placement. In the Voigt case, errors 
are dominated by continuum placement errors, and an estimate of this error is obtained 
by assuming a 1-a error in continuum placement, and then integrating this value over the 
wavelength range where the line profile is at a nontrivial depth (> 1%, a more than sufficient 
criterion for the deep lines under consideration). These two methods of error estimation are 
observed to be comparable when both are used on test cases. 

In most lines of sight, there is little to no evidence that the absorption profiles deviate 
significantly from a single component absorption structure, though the moderate resolution 
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of FUSE (approximately R = 20, 000) masks the true component structure^] However, 
in five lines of sight (HD 40893, HD 149404, HD 195965, HD 1999579, and HD 210839), 
there are consistently asymmetric profiles, which demand that we use a fit with multiple 
components. In these cases, we extend the methods described above to multiple components 
and record the total equivalent width. With this method we obtain an independent equivalent 
width for each component, but the equivalent width we use in the curve of growth is the 
total resulting equivalent width of the entire profile, not a linear sum of the two equivalent 
widths. We never find a consistent need for more than two components in these lines of sight 
that require multiple components. In the five lines of sight where multiple components are 
identified, we use the same curve-of-growth methods described below ( §2.5p . which employ 
varying assumptions. Table [3] lists all the equivalent widths and 1-a errors. 

The absorption lines of H2 are occasionally blended with other lines, bo t h oth er H2 



lines and atomic lines. The comparison list of atomic lines is from iMortoru ( l2003l ): HD 
lines are also considered. In these cases, our minimum criterion for proceeding with a line 
measurement is if the line centers are clearly resolved. In these cases, we simultaneously 
fit the partially blended lines and record the equivalent width fits. In cases where the line 
centers are unresolved or multiple components significantly complicate the fitting, we exclude 
these transitions from our analysis. 

In addition to obvious blending from stronger atomic lines or other observed H2 lines, 
there are many cases where the overlap between two H2 lines is nearly complete, such that it 
is not visually obvious that the feature is the blend of multiple lines. However, if the feature 
is clearly dominated by one line (e.g. a strong J = 2 line with a J = 7 line in a line of 
sight where J = 7 is otherwise undetected), it still may be appropriate for inclusion in our 
study. We determine the suitability of such lines by examining unblended lines (if they are 
detected) of the same J level of the line which is presumed to be weaker. Generally, if the 
presumed weaker line is from a state of J < 6, the stronger line is unsuitable for inclusion. If 
J > 7, then the suitability of the stronger line is dependent on that particular line of sight. 

A sample of absorption lines in HD 38087 is shown in Fig. [TJ This line of sight shows 
the greatest H2 excitation in our sample, as we detect up to J = 8. 



•^Note that a typical thermal width of H2 should be ^ 0.8 km s^^ and ^ 1.5 km s^^ for lower- and higher- 
J H2, respectively, based on the temperatures we derive in ^A.'2[ This calculation assumes a single cloud 
with no turbulence. These 6- values are significantly smaller than the FUSE resolution, which translates to 
a 6- value of ^ 9 km s~^. However, turbulent broadening or a multiple cloud structure may result in an 
observed 6-value that is comparable to the resolution. 



- 8 - 



2.4. Column Density Limits and Vibrationally Excited H2 

We also briefly explored the limits of the next undetected rotationally excited J level 
and searched for vibrationally excited H2. Our upper limits are based on the equation 



X,max S/N 

where Wx^max is the upper limit on Wx, is the number of a confidence desired, d\ is the 
wavelength spacing of the pixels, M is the number of pixels required to scan Wx, and S/N 
is the signal-to-noise of the local continuum. We assumed 15 pixels for detection (the FUSE 
FWHM is 9 pixels corresponding to a 6- value of ~ 9 km s~^), and calculated 1-cr upper limits 
on Wx'i column density limits can then be calculated from matching the limit on Wx to our 
preferred curve of growth solutions (see §3. 1.41) . 

While investigating these limits, we also detected a few lines of vibrationally excited H2 
in HD 38087 and HD 199579. The strengths and velocity offsets of these lines are consistent 
with our observations of the other H2 excitation levels. Detections of vibrationally excited H2 
in cold diffuse clouds are very rare, although at least two lines of sight from the FTCS sample 
that are not included in this study s how H2 in several o f these levels — HD 37903 (discussed. 



along with a review of the topic, in lMeyer et al.ll200ll ) and HD 164740 (B. L. Rachford et 



al. 2010, in prep.). However, because we only detected a total of five absorption lines in 
two lines of sight in this sample, we do not include these detections in the rest of our main 
analysis (Q. 



2.5. Derivation of High- J H2 Column Densities 

Few, if any, of the high-J absorption lines that we examined in this study are of an 
appropriate strength^ that a reliable column density can be determined directly from a 
single line. Rather, nearly all of the lines have at least some inherent saturation, but are not 
so strong as to show the clearly damped profiles of the J = — 1 lines (or allow for column 
densities to be accurately determined through the Wx oc \/N relationship that exists in this 
regime). Therefore, we use a curve-of-growth method to derive column densities. We use 
three versions of the curve-of-growth method in an attempt to be complete. 



''By "strength" we essentially mean the equivalent width, W\, which results from a combination of the 
line of sight properties (column density and velocity structure) and the inherent properties of the transition 
(oscillator strength / and damping constant 7). 
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The first method is to fit the fines of each J level of H2 to separate single-velocity 
dispersion curves of growth. Each J level is allowed to vary in column density and 6-value. 
The column density range is log = 10 — 22 in increments of 0.01 dex and the 6- value range 
is 6 = 0.1 — 30 km s~^ in increments of 0.1 km s~^, resulting in a 1201 x 300 grid of possible 
solutions in (iV, h)% This method avoids some assumptions that may not be justified, 
e.g. that different J-levels have the same 6-value, or that the velocity structure precisely 
matches the structure of some other atom or molecule. However, the uncertainties inherent to 
this method can be very large, both the statistical uncertainties when most of the absorption 
lines are on the "fiat" portion of the curve of growth, and the systematic uncertainties if the 
velocity structure is not closely approximated by a single- velocity dispersion. As a result, 
when we determine the error ellipse by calculating at each point in the column density/6- 
value grid, the resulting error is sometimes very large as well. In addition, we sometimes 
see unusual solutions (e.g. population inversions that are extremely unlikely to be real) that 
imply systematic errors. 

The second method is similar to the method described above, but instead we fit all the 
J levels of H2 to the same single- velocity dispersion curve of growth, instead of allowing for 
independent 6-values for each J level. This method helps to constrain the column density 
solutions, but whether or not this assumption accurately describes the physical conditions 
is uncertain. At the least, it prevents us from exploring one of the questions mentioned in 
^ concerning whether or not there is a difference in the 6- value of low- vs. high- J H2. 

The third method uses multiple-velocity dispersion curves of growth based on the veloc- 
ity structure of appropriate molecules. Specifically, we use the column densities and velocity 
structure of CH and CH+ that have been measured by D. E. Welty (2008, private commu- 
nication). This method is physically justified in that H2 and CH are known to be strongly 
correl ated (see the FTCS); fur thermore, higher-J H2 and CH"*" are both know n to be corre- 



lated (ILambert fc Dankslll986l ). We note that radiative transfer models (e.g. iNehme et al. 
20081 ) are unable to explain observations of the high- J H2 and CH+. Like the previous 
method of assuming a uniform 6-value, this method prevents us from directly examining 
whether or not low- and high- J deviate from full physical coincidence. Rather, this method 
makes an implicit assumption about whether this is the case, based on which J levels are fit 
to which species and whether CH and CH"^ are physically coincident or distinct. 

For all of these methods, the choice of a damping constant (7) affects the fit — this is 
virtually always the case for J = 2, occasionally the case for J = 3, and virtually never the 



^It is worth noting that the ranges on N and b are much broader than the expected sohitions, by > 2 dex 
for each hmit on N , and a factor of > 2 for the upper hmit on h. 
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case for J > 4. Therefore, in our calculation, the difference between the measured data 
points and the theoretical curves are based on different curves, appropriate to the damping 
constant of the absorption line under consideration. For example, the J = 2 line at 1081.267 
A has a damping constant of 71081.267 = 1-63 x 10^ s~\ while the 1016.458 A line has a 
damping constant of 71016.458 = 1-236 x 10^ s^^. Therefore, when calculating the value of 
for a column density of logiV = 18.50 and 6 = 2.0 km s~^, the difference terms in the 
sum use curves of growth with different values of 7, but the same 6-value. 

In order to simplify the computational process slightly, not every possible damping con- 
stant was calculated. The damping constants of the lines under considerations are bounded 
by 5 X 10^ s~^ < 7 < 2 X 10^ s~^; therefore, theoretical curves of growth were calculated for 7 
on this range, in increments of 5 x 10'' s~^. In the worst case, these values of 7 are separated 
by 10%, meaning that each absorption line will be matched to a value of 7 within 5% of its 
true value of 7. In the worst case scenario where lines are fully damped, W\Oi ^7; there- 
fore, this method is equivalent to introducing as much as a 2.5% error in VFau However, we 
note the following three mitigating factors: (1) the errors should be more or less randomly 
distributed over the dozens of lines measured, resulting in a reduced cumulative effect, (2) 
the introduced error is usually much smaller than this (it is negligible for the weaker lines), 
and (3) even in this worst case, this error is comparable to the statistical fitting errors on 
the measured values of and is usually smaller. 



3. RESULTS 

In this section, we will discuss the results of the first level of analysis based on the line 
measurements. First, we compare the various methods of deriving column densities. Then 
we will turn our attention to understanding the 6- values and velocity offsets of the absorption 
lines, and what this says about the nature of where higher- J II2 is formed and/or exists. In 
the following section we will use our adopted column densities to explore the additional 
issues of correlations with other atoms and molecular, excitation temperature, and radiative 
transfer modeling these lines of sight. 



^It should be noted, conversely, that W\ cx \fN on this regime, so the error in N would be as high as 5% 
if based on a single line. 
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3.1. Column Densities 

As discussed above, we performed column density analyses in three different ways — 
fitting column density and 6-value independently for each J level, fitting all J levels to a 
single 6- value (and independent column densities), and fitting all J levels to a predetermined 
curve of growth. Making quantitative comparisons between these methods is difficult because 
of the systematic uncertainties. Nevertheless, we must attempt to understand which of these 
methods produces the most accurate results. 

3.1.1. Independent b-values 

When we assume that 6- values are independent, we notice, in general, trends of in- 
creasing b with increasing J. Compared to a uniform 6-value, this will tend to decrease the 
column densities of the higher- J H2. It is important to note that in some but not all cases 
this effect is neghgible because the error in b is large — the best fit is simply finding that the 
higher- J H2 appears to be on the linear portion of the curve of growth, and the 6- value fit 
and its error should instead be interpreted as a lower limit. 

We have attempted to look for direct causes of potential systematic uncertainties by 
examining anomalous absorption lines. There are certain lines which frequently have equiv- 
alent widths that are much larger or smaller than expected; lines which are systematically 
stronger or weaker than expected across multiple lines of sight are consistently identified 
as having a continuum that is less certain than other lines. When we exclude these lines 
that have apparent systematic uncertainties from the curve-of-growth analysis, the curve-of- 
growth fits are improved in many cases. However, even after this consideration there are still 
some column density results that we consider anomalous — specifically, strong deviations from 
a reasonable population expectation, that is, population inversions isolated to one J level. 
Possible additional systematic uncertainties that we have not considered include undetected 
blending with other atomic or molecular transitions or /-value errors. 

Our conclusion is that this method, while preferable aesthetically because of its lack 
of a priori assumptions, is ultimately the least reliable, because there are fewer constraints, 
namely a much smaller range in /A. This enhances the potential for systematic errors in the 
final column density derivation due to inherent deviations from a single-component curve of 
growth or unrealized systematic errors in the measurements of the absorption lines. 



- 12 - 



3.1.2. Uniform h-value 



Assuming a uniform 6-value for all J levels of H2 is a good way to put constraints on 
the curve of growth method and obtain a seemingly reasonable solution. However, while a 
uniform fe-value seems like a good first-order assumption, there is a distinct possibility that 
this may not b e the case. Discussio n of this p henomenon can be fo und in earlier Copernicus 



result s, such as lSpitzer et al.l (Il974j ). Papers by I Jenkins et al.l (Il989l ) and I Jenkins &: Peimbert 



(119971 ) also found some evidence that 6- value varies as a function of J — although the nature of 
the effect was the opposite for the two lines of sight in those papers (decreasing 6-value with 
increasing J the former paper's examination of vr Sco, increasing fe-value with increasing 
J in the latter pap er's discussion of C Ori A). A more recent example can be found in 
Lacour et al.l (120051 ) . who determined to 4-cr confidence that HD 192639 cannot be fit by a 
single 6-value — h increasing with increasing J in this case. In their paper, iLacour et al.l also 
discuss the implications this has for the formation mechanisms of H2. This line of sight is 
also in our sample, so we are able to perform the analysis using both methods (independent 
and uniform 6-value) on this line of sight. When fitting 6- values inde pendently to e ach J 
level, we do not rule out a single 6-value at the same high confidence as iLacour et all J2005h 
but our results are qualitativ ely similar. In addi tion, quantitative differences between our 
column densities and those of ILacour et al.l (120051 ) are observed at levels between 1- and 3-0" 
for J = 2 - 4. 

However, it may still be useful to attempt to fit all J levels to a single fe-value. As a 
practical matter, this helps to restrict 6-value solutions to qualitatively reasonable values 
(i.e. no population inversions), but may introduce systematic errors if this assumption is 
invalid. Furt h ermo re, in spite of the aforementioned potential evidence to the contrary in 
Lacour et al.l (120051 ) . assuming a single 6-value is the simplest a priori assumption, and is 
worth more exploration. 

Even within the context of assuming a "uniform" 6-value, we might still find multiple 
6-value regimes as a function of J. What if, for example, H2 in the J < 3 states is at 
one 6-value, while H2 in the J > 4 states is at another? This would be consistent with 
the observation of different excitation temperature regimes as a function of J (see §4.2p . 
Unfortunately, we cannot probe the 6-values of J = — 1 because these lines are strongly 
damped and therefore the profiles are less sensitive to changes in 6-value. Therefore, between 
continuum uncertainties and the ~ 15 km s~^ resolution of FUSE, profile fitting cannot easily 
differentiate between 6-values for these lines. Curve-of-growth methods likewise cannot easily 
distinguish between 6-values, with all the lines being well onto the damped portion of a typical 
curve. This is also potentially a problem for the J = 2 lines, which largely exist at the flat- 
to-damped transition and are only mildly sensitive to 6-value — however, with the stronger 
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J — 2 lines setting the column density, the weaker J — 2 lines are frequently able to help 
set the 6- value with reasonably high confidence. 

To explore this, we performed three variations of the "uniform 6-value" method: 

1. Assume that b is uniform for J > 2 (62+) 

2. Assume that b is uniform for J > 3 but independent for J = 2 (63+ and 62, respectively) 

3. Assume that b is uniform for J > 4 and an independent 6- value is shared by J = 2 — 3 
(64+ and 623, respectively) 

First, we note that there is a great deal of consistency for J > 4 for all three methods — 
16 of the 22 lines of sight match within their 1-a column density errors. An additional 
three fines of sight (HD 27778, HD 179406, and HD 206267) match within 1-a errors for 
J > 5 column densities, but have as much as a 0.04 dex gap between the largest lower limit 
and the smallest upper limit for J = 4. The remaining lines of sight are HD 24534, HD 
53367, and HD 210839. The discrepancies for HD 210839 are larger than the errors but not 
extreme: 0.16 dex for J = 4 and 0.06 dex for J = 5. For HD 24534 and HD 53367, aU of the 
discrepancies are at least 0.24 dex, and are greater than 1 dex for the two HD 24534 column 
densities. In these cases, methods 1 and 2 are consistent, but method 3 sets values of 64+ 
and 623 that are significantly different both from each other and from the 6-values derived 
by methods 1 and 2. 

In general, the biggest variations are for J = 3. While several of our lines of sight are 
either consistent within their errors or have minimal variation (5 agree within errors; another 
3 agree within errors plus 0.04 dex), the majority of lines of sight are not this close, with 
disagreements of over 1 dex in many cases. In 20 of 22 lines of sight, method 3 is the outlier 
of the three methods, while methods 1 and 2 are consistent within errors for 15 lines of sight, 
and consistent within errors plus 0.1 dex for another five lines of sight. 

The nature of these discrepancies makes sense because the J = 2 equivalent widths, 
nearing the damped portion of the curve of growth, are less sensitive to fe- value and therefore 
do not constrain it as much in the method 3 fit; the J — 3 lines are then free to determine 
the 6- value in the fit, and suffer from the systematic uncertainties similar to the independent 
6- value case. However, in methods 1 and 2, where the J — 3 equivalent widths are "linked" 
to the J > 4 equivalent widths, the J — 3 column density fits are constrained to a better- 
justified 6- value. These qualities cause us to rule out method 3 compared to the other two 
methods. 

There is also some variation in the J — 2 fits between the three methods, greater than 



- 14 - 



the variations in J = 4 — 5 but not as much as for J = 3; again, this makes sense because 
the J = 2 column density fits should be less sensitive to 6-value. Comparing methods 1 and 
2, where J = 2 is linked with the other column densities versus allowed to be fit to its own 
6- value, we see consistency within the errors in the majority of cases (12 of 22). Another 
six cases have disagreements of 0.07 dex or less (typically 2- or 3-a discrepancies). In the 
remaining cases, a qualitative check favors method 1 — adopting method 2 would imply a very 
small ratio of J = 2 to J = 3, to the point of having a negative excitation temperature, T23 
in two of the four cases (see also §4.21 for more general discussion of excitation temperatures). 

The consistency check between these variations gives a sense of what potential system- 
atic errors are inherent to the different methods. We have concluded that method 1 is the 
best version of these three variations, based on the simplicity of its underlying assumption 
and the avoidance of unusual solutions (e.g. a negative T23). We note, however, that in 
principle any individual line of sight might be better fit by a different variation. Further 
discussion concerning the best choice of curve-of-growth methods is found in §3.1.41 



3.1.3. Predetermined Curves of Growth 

We calculate curves of growth from fits of CH and CH+ column densities made by 
D. E. Welty (2008, private communication) using various high-resolution (1.2 — 3.6 km s^^) 
ground-based optical data set^, and then fit the H2 lines to these curves. In creating these 
curves of growth, the 6-values have been scaled with the mass of the molecule, according to 
the formula: 

h^ = '^-^+e (2) 

m 

In this equation, m is the mass of the molecule and ^ is the turbulent velocity. Because 
the turbulent velocity and temperature are not explicitly known, one or the other must be 
assumed. Therefore for each identified cloud component of CH or CH"*" with a measured 
6-value, we assume a temperature, then use the 6-value and the molecular mass to solve for 
^. We then use that value of ^ and the H2 mass to calculate the new assumed 6-value for H2 
in that cloud. 

The temperature that is assumed is 70 K for all lines of sight, which is broadly consistent 
with the average core temperature that was derived in the FTCS and our study (see §4.2p 
through measurement of the low-J excitation. However, the variations in 6-value due to 
inaccuracies in this temperature assumption are small — h varies as ~ a/T for small ^, and as ^ 



^Information about the majority of the data can be found at 



http: / / astro.uchicago.edu/ ~welty/ 
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becomes larger, small differences in temperature have a reduced effect. In our data, ^ usually 
dominates the 6- value, i.e., ^ is usually greater than \J^^-- Therefore, universally assuming 
a temperature of 70 K should result in a reasonable scaling of the 6-value. Attempting to 
further model the temperature will not result in significantly increased precision due to the 
inherent uncertainties of the assumptions that are being made. Another issue is that the 
higher J levels are found at a higher excitation temperature (again see 14.21) . but this may 
or may not translate to the thermal temperature because the higher J may be partially 
influenced by radiation. With the various combinations of CH, CH"*", and lower- and higher- 
J H2, we would potentially have to assume multiple temperatures with little observational 
guidance. We also note again that when ^ > \J^^i even a factor of a few in temperature 
may not have a significant impact on the scaled 6-value. However, scaling the 6-value in 
this manner with the assumption of a 70 K temperature should provide an improvement 
compared to no scaling. 

The foremost question is whether to use the CH or the CH+ data. As discussed above, 
total H2 column density [dominated by A^(0) and A^(l)] is generally correlated with CH and 
higher- J H2 is correlated with CH"*". However, it would be prudent for us to confirm these 
correlations in our data set, and also note that correlations do not necessarily mean that 
assuming a common velocity structure is the most appropriate curve-of-growth method. 

Attempting to both look for correlations and assume a correlation so that a curve of 
growth can be assumed involves obvious circular reasoning. Below in §4. II we have examined 
the different J-levels for correlations with CH, CH"*", and other atomic and molecular species. 
We find that some of the most consistent correlations are with CH and CH+. This was done 
for the H2 column densities determined assuming a uniform 6-value. We do confirm a high 
level of correlation for both CH at very low- J and CH"*" with J = 3 — 5. However, neither is 
correlated at statistically significant levels for J = 2 or J > 6. It is also worth noting that we 
only have velocity structure data for 17 of our 22 lines of sight, and both CH and CH+ data 
for only 14 lines of sight (note that there are additional column density derivations, which 
are based on equivalent width measurements but not a full fit to the component structure, 
that are available and included in our analysis in §4.11) . Therefore, if we were to adopt this 
method, we would not be able to perform a uniform data analysis over the whole set. 



3.1.4. Column Density Conclusions 

In Fig. [2] we show an example of how variable the different methods are. This is done for 
HD 38087. In this line of sight, A^(2) and A^(3) are completely consistent within their errors 
[note that N{0) and A^(l) are shown in the plot but not independently derived]. From A^(4) 
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through A^(6), and to a lesser degree N{7), there is some disagreement between the methods; 
at A^(8) the methods are once again consistent (noting the very large error on the independent 
6- value derivation for N{8), because the limited number of absorption lines are consistent 
with a wide range of 6- values using the curve of growth method). The inconsistencies are 
obviously due to the differing levels of saturation that is implied by each method. What 
may be initially perplexing is that the column densities based on the CH velocity structure 
measured in two different data sets have the most variation (again when the large errors of the 
independent 6- value results are considered). This is explained by the inconsistency between 
the 2000 and 1991 data sets for CH. The fit 6-value of the single observed CH component in 
these two data sets differs by a factor of nearly two; however, the fitting errors in the 2000 
data set are much smaller (2.2 ±0.4 km s~^ compared to 1.3 ± 1.8 km s^^), because the 1991 
data set is particularly noisy (D. E. Welty 2009, private communication). The much smaller 
6-value implies that the moderate-strength absorption lines of the J = 4 — 7 are much more 
saturated, hence the larger resulting column densities. Because our fixed curves of growth 
are calculated without consideration of the error in the fits in these predetermined 6-values, 
we cannot account for uncertainties such as this in any way that is not computationally 
and conceptually complex. It is somewhat notable how consistent the independent 6-value 
method is with the variations on the uniform 6-value method within their errors, with the 
exception of A^(6). 

Overall, we find a fair amount of consistency between the three methods. The more con- 
sistent methods are assuming a uniform 6- value and using predetermined curves of growth — a 
signal that using a uniform 6-value is frequently a fair approximation of the true velocity 
structure. However, there are occasionally significant discrepancies both between the three 
methods (independent 6- value, uniform 6- value, and fixed velocity structure), and within 
their variations (assuming a "break" in the 6-value of the uniform b method or different fixed 
curve-of-growth component structures). 

We have elected to use the results from the uniform 6-value method, assuming that the 
6-value is constant across all J levels. This is a decision that is equal parts quantitative 
and qualitative. Despite the possibility that 6-values may vary as a function of J (see 
earlier discussion in §3.1.21] . we cannot confirm this across multiple lines of sight (see §3.2.ip . 
Furthermore, while we find CH and CH"*" to be strongly correlated with different J levels of 
H2 (see §4.11) . the correlation is weak for J = 2 and J = 6 — 7. Even the strength of the 
correlations for the J levels which are correlated do not guarantee that the velocity structures 
are close enough to ensure this method is the most accurate. Finally, we note that our data 
set would not be uniform in this case, because we do not have complete CH and CH+ data 
for all of our lines of sight. 
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Quantitatively, the major argument against using the measured velocity information 
from CH and/or CH+ is that it rarely improves the compared to the uniform 6- value 
method, never does so by more than 10-20%, and in several cases results in a significantly 
worse (factors of a few, up to ~ 7 in the case of HD 199579). Similarly, the quantitative 
argument in favor of the uniform 6-value method is that the value it produces is never 
worse than ~30% of the best x^ of all methods we explore, with differences of only a few 
percent being typical (and it is the best method for some lines of sight). Furthermore, because 
of its nature, we do not see any unlikely population inversions or dramatically inconsistent 
6- values. 

Figs. [3111] show the curve of growth fits for all 22 of our lines of sight assuming a uniform 
6- value. The fits to the various lines are color-coded by J level. Table S] provides the column 
density results. Note that the errors reported in this table, for computational simplicity, 
assume that the other column densities and 6-value are held fixed, though technically the 
fitting of these column densities is dependent on the other column densities and their influ- 
ence on the 6-value. Limits of the first undetected J levels and vibrationally excited levels, 
and the few detections of vibrationally excited material, are shown in Tables [MEl 

Though we adopt the uniform 6- value method, comparison of the CH and CH"*" velocity 
structures with our adopted results are still instructive (we compare column densities in 
§4.ip . For example, the five lines of sight where we derive the smallest 6- value and also have 
CH and/or CH+ velocity structures (HD 24534, HD 27778, HD 38087, HD 147888, and HD 
170740) all have 6-values of < 3.1 km s~^. These same lines of sight have only one or two 
detected CH/CH+ components, separated by < 3 km s~^. Only HD 210839, where we solve 
for a 6- value of 9.7 km s~^, has a comparably small velocity extent in CH"*" (~ 2 km s~^), 
but the CH, K I, and Na I velocity structures show many more components with wider 
spacing in this line of sight. Conversely, the largest 6-value that we solve for, 12.5 km s^^ 
in HD 199579, matches with the CH velocity profile with the widest spread of components 
(~ 22 km s"^); HD 149404, the third highest 6-value (after HD 199579 and HD 210839) for 
lines of sight with CH/CH+ data, is similar. The remaining lines of sight with intermediate 
values for this study (~ 4 — 8 km s~^) are likewise intermediate in terms of the complexity 
and velocity extent of their CH/CH+ profiles. 

However, note that this is a correlation of the velocity structures and not a strict corre- 
spondence. In fact, as the structures become more complex, some of the fits to the observed 
CH/CH"*" structure are very poor. Fig. [5] shows the case of HD 210839 fit to the CH velocity 
structure. In this case, there is a huge gap between the derived column densities of J = 5 and 
J = 6, and the fits for these column density levels are particularly poor. Compare this to the 
last panel of Fig. HJ where, while there is some scatter, the points follow a something close 
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to single-component curve of growth (even though equivalent widths are of the measured 
profiles that appear to have at least components even in the FUSE data). The situation is 
similar for the CH"*" velocity structure in this line of sight, and the CH/CH+ structures in 
some other lines of sight, particularly HD 199579. 

3.2. Velocity Structure as a Function of J 

As we noted in ^ the excitation of H2 is so poorly understood that certain first-order 
types of questions about its velocity structure need to be answered. For example, it is possible 
to construct a qualitatively simple model of a line of sight where, due to the excitation 
processes of H2 and the physical structure of the cloud(s), different J levels of H2 essentially 
exhibit a distinct velocity structure. For example, there might be two major clouds in a line 
of sight, one warm and one cold. In this case, the former cloud might dominate the observed 
lines of higher- J H2 and the latter cloud the lower- J lines. Another simple possibility is that 
the higher- J H2 exists in a warm shell around the cold core containing the low- J H2. In both 
of these examples, a possible (but not necessarily inevitable) observational consequence is 
that the simplest measures of velocity structure — 6- value and velocity offset — may vary as a 
function of J. Therefore, we conducted an examination of both of these quantities. 

3.2.1. b-values 

There are two basic ways to examine whether or not 6-values vary as a function of J. 
The first is to look at the apparent 6- values of the Gaussian and Voigt fits. The second is to 
examine the 6-values derived through the curve of growth fits. 

In examining the apparent 6- values of the fits, we use only the 17 lines of sight for which 
profiles were fit with a single component. We use a simple linear regression with Ej/k as the 
abcissa, average 6- value for a given J level as the ordinate, and the standard deviation in the 
6- value for that J level as the error. We do this for each data channel of FUSE separately to 
provide independent checks on instrumental effects. The slopes are nearly always negative 
and typically significant at between 1- to 3-a (there are two isolated exceptions of positive 
slopes in single data segments that are less than l-cx significant; in both cases, the other 
three data segments have negative slopes). There are also a few cases of slopes that are 
more significant than 3-a. We calculate correlations with IDL's R_CORRELATE function, 
which computes a Spearman's p rank correlation coefficient. A rank correlation is used in 
this case because we do not expect a particular functional form for any dependence. Rather, 
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we want to search to see if a roughly monotonic relationship exists — e.g. 6-value increasing 
as a function of J . We also compute the probability of deviation from the null hypothesis, 
which gives similar results — nearly every data segment shows that an anticorrelation exists, 
typically significant at the 1- or 2-cr level. 

However, we should note that when we analyze the unresolved apparent 6-values of the 
fits, they are obscured by the broad resolution of FUSE (~ 15 km s~^). This complicates 
the matter in two ways. First, when the fits are Gaussian, the apparent width increases as 
a function of column density due to natural broadening, while still appearing very nearly 
Gaussian in the FUSE data. Second, if there are weak components separated in velocity 
space, as these components become stronger with increasing column density the apparent 
6-value will also increase — that is, these components are lost in the noise for weak lines but 
blend with the domi nant component (s) for strong lines, an effect that was noted in early 



studies of the ISM by lRoutly fc Spitzerl ( 1l952l ). Because the lower J states have larger column 
densities, both of these effects could result in the negative slopes that we see — which say 
that 6-value is decreasing as a function of Ej/k. In addition, the Ff/S'i? resolution varies as a 
function of wavelength within individual detector segments, which introduces an additional 
level of uncertainty. Together these effects significantly reduce our abihty to use this method 
to see any real cases of 6-values changing as a function of J, because at the resolution 
we cannot adequately account for these effects. Fig. |6] shows histograms of both correlation 
coefficients and linear regressions, showing the systematic offset due to these effects; it can 
be seen in these histograms that additional deviations are likely no more than statistical 
variations. 

Looking at the 6-values derived from the curve of growth fits suffers from similar prob- 
lems. The apparent width due to natural broadening is not an issue (because the width is not 
being directly measured) , but the possibility of an apparent increase in 6- value with increas- 
ing cohjinn_d^n^ity_(ai^^ lower J value) due to the weak component effect noted 



by iRoutly &: Spitzerl ( 119521 ) still is. In addition, the systematic uncertainties of assuming a 



single- component curve of growth come into play. 

Again using IDL's R_CORRELATE function, we find that 12 of 22 fines of sight show 
correlations are less than l-a significant; the other 10 are less than 2-cr significant. Exactly 
half of each level are correlations; the others are anticorrelations. From this, we might con- 
clude that any deviations are statistically distributed fiuctuations. When we use linear fits 
(6- value vs. Ej/k), there are several (8) lines of sight that have greater than 3-a signifi- 
cant slopes. However, these are mixed both positively and negatively. The reason for the 
high confidence in the slopes may come from an underestimation in errors in the 6-values. 
The mixed signs of the slopes and correlations mean that we do not see a clear universal 
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dependence of 6-value on J (in particular, we might expect an increase of b with J due to 
formation processes). This is particularly true given the potential systematic errors that we 
have discussed, though at the same time weakens the confidence with which we can state 
the null result. Fig. [7] shows histograms of the correlation and regression methods, showing 
that the deviations may be explained as roughly statistical variations. 



3.2.2. Velocity Offsets 

We also examined velocity offsets as a function of J in each line of sight. We take an 
average of the velocity offsets of all the lines of a given J level in each line of sight and in 
each FUSE channel. Keeping the FUSE channels separate is more important for velocity 
offsets than for 6-values because there are systematic offsets between the channels which 
our initial data reduction does not correct (it is relatively unimportant other than for this 
question). Again, it also allows us to perform a check on the significance of any possible 
trends — a slight trend may be seen in one data channel but not in the others. As in §3.2. 
we only perform this analysis for the 17 lines of sight where we used single-component fits 
to measure the velocity offset. 

One of us (BLR) measured velocity offsets for the J = — 1 lines from previous data as 
part of the FTCS, but there are problems with trying to use these data in our comparisons. 
First, some of these fits were made with very early versions of the CALFUSE pipeline reduc- 
tion that may not be consistent with ours, which has the potential to be very significant for 
measuring velocity offsets. Second, these previous fits included the J = 2 lines, meaning that 
the velocity measurements of the J = 0, 1 lines are not completely independent. Therefore, 
we instead measure the HD lines that are in the FUSE wavelength range. This is in part 
justified by the very strong correlation that we find betwee n HD and H 2 (see jj4.1l) . We use 



the same HD lines discussed in the recent survey of HD by ISnow et al.l (120081 ) . except that 
we omit the weakest line at 1105.834 A because it is generally too weak to detect and when 
it is detected it suffers from confusion with neutral carbon lines. Because we measure the 
J = lines of HD, we use these as a proxy for H2 in the J = state and omit J = 1 from 
this portion of the analysis. 

When we examine only J > 2, there are few correlations to be seen. There are no cases 
where as many as three of the four data segments for a given line of sight have greater than 
1-cr significance (judging by slopes); occasionally 1- or 2-cr significant correlations appear 
in one or two individual segments, but are not supported by the remaining two or three 
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data segments. Plotting a histogram of the 68 correlation levels and 68 slope significanceqj 
indicates that any exceptions are easily explained as statistical variations. 

However, when the HD lines are included as a proxy for J = 0, the situation becomes 
more complicated. The rank correlation coefficients are centered around a positive l-cr 
deviation, and the slope significances are centered around negative 1.5-(j|^ Though this 
superficially implies that there is a discrepancy between the location of the coldest H2 (where 
the HD is) and the higher- J H2, there are issues to consider. First, if there is such an 
effect, we would not expect it to be manifested as a shift toward either positive or negative 
correlations that holds across all lines of sight, we would instead expect that there are 
more outliers that can be explained statistically when assuming a null hypothesis. The 
width of the correlation/slope significance histograms stays constant in our case. Next, 
we measure a relative paucity of HD lines compared to H2 lines — we only search for six 
lines, and in some cases only one or two lines are measured due to poor S/N. Therefore, 
it is possible that there are systematic errors in the measured velocities of the HD lines, 
either in our reference wavelengths or due to the instrument. We must consider that the 
FUSE resolution is 15 km s"^ (FWHM) and that there may be comparable uncertainties 
in the FUSE wavelength solution. While our measurements of many lines are intended to 
reduce these errors through improved statistics, they are still significant. The average slope 
measured in the HD-inclusive case, integrated between the J = HD and the maximum 
J-level of H2 detected for each line of sight, is 10 km s~^. Finally, extending the correlation 
down to J = 2 should be sufficient to detect differences between the coldest and warmer 
material, as A^(2) is consistent to a degree with the excitation temperature implied by iV(0) 
and iV(l), further supporting the possibility that the HD lines introduce systematic errors. 
Fig. [8] shows histograms of the correlation coefficients and regression significances. 

Ultimately we conclude that we do not observe any definite cases of velocity offset 
varying with J, as the HD-inclusive velocity results, while intriguing, are more likely due to 
systematic errors. It is very important to note that the imprecision of the Ff/^i? wavelength 
solution and its moderate resolution prevent us from stating a null result more strongly. 
However, our results are still suggestive of the conclusion that the higher- and lower- J H2 
are not dominated (in an observational sense) by physically distinct regions. 



^Four data segments for each of the 17 lines of sight that use only single-component fits to the absorption 
lines. 

^The correlation coefficients are translated into probability deviations; positive deviations correspond to 
negative correlations/slopes. 
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4. DISCUSSION 

We now turn our attention to other important issues regarding the interpretation of our 
measurements: how does higher- J H2 correlate with other molecules, what is the temperature 
of these clouds (and what is the implication for formation and destruction processes), and 
can these results be modeled theoretically? 



4.1. Correlations with Other Atomic and Molecular Species 

CH is known to be correlated with total H2 (dominated by J = — 1; see the FTCS 
and references th erein) . In lines of sight with less extinction, CH"'' is known to correlate 



with higher- J H2 (iLambert fc Dankslll986l ). What has not been explored in detail at higher 
extinction is how higher- J correlates with CH, CH+, or other atoms and molecules. We 
have taken available column densities for several of these atomic and molecular species and 
attempted to find existing correlations. 

Our original working list of atoms and molecule^ included C I, S I, K I, Li I, Ca I, 
Ca II, Na I, Fe I, Fe II, Ti HI, C2, C3, ^^CO, ^^CO, CN, HD, CH, and CH+. Obviously, 
only some of these are expected to be highly correlated with H2 regions, while others may 
give an indication of the radiation levels. Notably, the database we are drawing from does 
not overlap perfectly with ours, meaning that in some cases we are attempting to correlate 
a very small number of points. Though we calculate the probably that correlations are not 
null based on the number of available data points, any correlation with a significantly small 
number of data points should be considered with healthy skepticism. We have therefore 
chosen to only include in our calculations those atoms and molecules for which there exist at 
least 10 lines of sight overlapping with our study (and therefore > 10 common data points 
for all correlations involving J < 5; correlations involving J > 6 may have fewer points). 
This revised list is Li I, K I, Ca I, Fe I, Fe II, HD, CH, CH+, ^^CO, C2, and CN. 

Tables WM show correlations between these species and various measures of the H2 
column density and population distribution. All the tables show the Pearson correlation 
coefficient r (positive is a correlation, negative is an anticorrelation); the two-sided confidence 



^"These data are primarily taken from a database of colmiin densities by D. E. We lty available at 

Many of these data are published in IWeltv fc Hobbs 



ht t p : / / www . astro . uiuc . edu / ~ dwelty / coldens . ht ml 
( 200lJ )^ Weltv et al ) feOOB*^ : ISonnentrucker et al. ( 2007 ). Exceptions are that, where available, Fe II is taken 



from [ Jensen fc Snow, ((2007il and HD is from the single-component curve of growth results from iRoss et al.l 
|2008^ 
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level that the correlation is nonzero for a t-distribution of the quantity r(r?, — 2)^/^(1 — r^)~^/^ 
for n — 2 degrees of freedom, where n is the number of data points; and finally, n is also 
listed for each correlation calculation. 

Tables [71|8] have the same columns but different rows. The first column is total H2 column 
density, and the rest are several different ratios of the column densities of different J levels. 
The various ratios probe different qualities of the lines of sight. A^(2)/A^(0) and A^(3)/A^(l) 
are indicators of the g as density; greater ratios imply lower density, because J = 2 — 3 are 



colhsionally deexcited (IShuU k BeckwithI Il982f l. Similarly, A^(4)/iV(0) and iV(5)/A^(l) are 



indicators of the radiation field, because J = 4 — 5 column d ensities are set largely b y the 



radiation field, so larger ratios imply more radiation (again, see lShuU &: Beckwithlll982[ ). The 
other three ratios that we include in this table are A^(4)/iV(2), iV(5)/iV(3), and A^(6)/iV(2). 
These measures in essence merge the more commonly used "radiation" and "density" ratios 
that we just discussed to give combined measures of the two quantities; the last quantity, 
A^(6)/A^(2) does so by looking at A^(6) to examine the top end of the radiation, but is limited 
by the number of J = 6 detections. Table [7] shows the correlations between the parameters 
discussed in the previous paragraph with the column densities of the atomic and molecular 
species. Table E] explores three potentially interesting ratios of the atomic and molecular 
species: CH+/CH, CH+/CN, and Fe II/Fe I. 

In Table [71 the correlations for Li I, Ca I, Fe I, Fe II, C2, and CN are less than < 95% 
significant with only the correlations between total H2 and C2 and CN > 90%. K I shows 
a significant correlation with total H2, and significant anticorrelations with the density and 
radiation ratios. All of these are unsurprising, as we expect K I, with its lower ionization 
potential, to only exist were it is relatively dense and shielded. HD shows a similar pattern, 
although the correlations for the density and radiation ratios do not meet the criterion for 
significance as they do for K I [although A^(5)/A^(l) is close at 94.7%]. 

CH, as already known, is very strongly correlated with the total H2, but it is almost 
as strongly anticorrelated with the combined ratio of A^(6)/A^(2). Anticorrelations are also 
observed for the density and radiation ratios, although only A^(3)/A^(l) meets our significance 
criterion. For CII+, significant correlations are observed with the radiation ratios, and two of 
the combined ratios are significant or nearly so [94.5% for A^(4)/A^(2), 97.7% for A^(5)/A^(l)]. 
Combined with the fact that the correlations with the density ratios are insignificant, this 
stresses the possible importance of radiation in CH"*"; at the very least, it confirms the 



^^That is, the probability that the calculated correlation coefficients are not from a normally distributed 
population is < 95%. In the discussion that follows, 95% is the nominal criterion for "significance," although 
some exceptions are discussed. 
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connection between CH+ and J = 4 — 5, including potentially between their formation 
processes. 

Finally, ^^CO has a correlation pattern that is basically the opposite of CH+ — significant 
anitcorrelations for A^(4)/A^(0) and A^(4)/A^(2) and nearly significant anitcorrelations for 
A'"(5)/A^(l) and A^(5)/A^(3) (just over 90% but not meeting our significance criterion). As 
photodissociation is the dominant destroyer of CO, these anticorrelations are not surprising 
to the extent that the radiation ratios and combined ratios do actually correspond radiation 
levels. 

Moving on to Table IH of the three ratios shown in the rows, only CH+/CH shows signif- 
icant results. The correlations with the radiation ratios are very strong (> 99.5%), as would 
be expected based on the correlations with CH and CH"*" separately. The other correlations 
also follow according to expectation, although at reduced significance — the combination ra- 
tio iV(4)/(2) is still significant at 99.0%, while the combination ratio A^(5)/A^(3) and the two 
density ratios are > 90% but do not meet our significance criterion. Once again, we note 
the strong connection to J = 4 — 5 and, implicitly, to the radiation levels. 

Table [9] shows column density correlations as a function of J-level. In this table, 
the correlation coefficient is modified to account for the total hydrogen column density 
[A^ (Htot) = A^(H I) + 2 A^(H2)] of these hues of sight. That is, we adopt the formalism used 



by iJenkins et al.l (Il986[ ) in order to make sure that correlations are not the spurious result of 
total column density — if this step is not taken, all of the correlations are very, very strong. 
We also note that in doing this, there may be some effect of inhomogeneity (and therefore 
possible systematic errors) of the atomic hydrogen sources. The references for all of the 
atomic hydrogen va lues that we use , along with some discussion of concerns of inhomogene- 



ity, can be found in lJensen fc Snowl ( 120071 ). 



Many of the correlation coefficients of Table [9] are statistically insignificant, but there are 
exceptions. All the J levels are anticorrelated with Fe I, with the confidence levels meeting 
our selection criterion for J = 2 — 5. This is most likely explained as an effect of radiation 
ionizing the Fe I to Fe II. K I follows the correlation patterns seen in Table [7] (correlated for 
J = — 1, anticorrelated for J > 2). Only J = 2 and J = 6 meet our significance criterion; 
however, all other levels except J = 7 are > 90%. 

The molecules are likewise consistent with what we would expect. HD is correlated with 
J = — 1 and anticorrelated with J = 5 — 6; only J = meets our significance criterion, 
but the other three are > 90%. CH is strongly correlated with J = — 1 (meeting our 
criterion) and weakly anticorrelated with J = 6 — 7 (failing to meet the criterion in part 
because of fewer data points; the J = 7 coefficient is actually quite large). CII+ is correlated 
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with J = 3 — 5 (at 94.9%, 98.9%, and 99.9%, respectively); none of the other correlations 
or anticorrelations meet the significance criterion. And lastly, the correlation coefficients 
roughly increase as a function of J for ^^CO — correlated for J = and anticorrelated for 
J > 4. Although only J = 5 meets the significance criterion, J = 0, J = 4, and J = 6 are 
all > 90% ( J = 6 is less significant than J = 5 despite having a larger correlation coefficient 
due to fewer data points). 



4.2. Excitation Temperature 

A well-established observation of H2 is that the column densities of the various J levels 



are r arely able to be fit by a single excitation temperature (e.g. see ISpitzer fc Cochran 



I973I ). alth ough exceptions are found in lower-column density, higher-temperature lines of 



sight (e.g. iMorton fc Dinersteinlll976l ). Excitation temperature is defined by 



Ni gi 

where N is the column density, g is the Gaunt factor (statistical weight), the u and / sub- 
scripts are the upper and lower state values, AE is the difference in the energy levels (upper 
minus lower), and Texc is the excitation temperature. If the logarithms of the statistically- 
weighted column densities is plotted against AE/k, then this relationship should be linear 
for constant Tg^c- However, typical plots show that this relationship is never linear from 
J = extending to J = 5 and beyond. Rather, most lines of sight are best fft with a sum 
of populations at different temperatures. Although the best fit may require more than two 
temperatures, it is typical to use one temperature from J = 0toJ = 2or3, and another for 
higher J-levels. Because of number of our data points (typically only 6 or 7 J-levels) and 
the number of free parameters required (one temperature and one normalization for each 
population assumed), we will limit our exploration to only two temperatures. 

Our first step is to determine the most appropriate break point for this fit. Defining 
Jbreak as the J Icvcl whcrc we assume J < Jbreak is fit with one temperature and J > Jbreak 
is fit with another, we test Jbreak = 2, 3, and 4. In every line of sight, Jbreak = 3 produces 
the lowest global reduced (= xl) — although xl is rarely below 1, revealing that either 
errors have been underestimated or this simple model is inadequate. 

Figs. [9lfT0] and Table [10] show our fits assuming Jbreak — 3. For comparison, Tqi-, derived 
solely from the J = — 1 column densities, is also shown. It can be seen that T02 and Tqi 
are not entirely consistent — only five lines of sight agree within their 1-a errors. However, 
the discrepancies are still instructive — in the five worst cases, the ones which disagree by 
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more than 3-cr, T02 > Tqi. This is what would be expected if the processes responsible for 
populating the J > 3 levels of H2 also have an impact on J = 2. 

We also note that for 19 of our 22 fits, the A^(3) is greater than the fit; in 14 of 
those cases, the discrepancy is greater than 1-a. Combined with what was discussed above 
about the difference in Tqi vs. T02, this is evidence that A^(2) and A^(3) are best fit by a 
third temperature. However, it is difficult to put quantitative certainty on this because the 
number of data points is too small, as discussed above. We also note that the relationship of 
A^(4) and N[5) to the fit do not seem to indicate an ortho-para effect; both N{4) and N{5) 
tend to be underestimated by the fit greater than half of the time, although several points 
are within the errors. 

Judging by the plots, it is also not infrequent that the column density of the highest 
J level commonly detected — J = 6 — is underestimated by the fits. There is certainly the 
possibility that H2 in the J > 6 states is better fit by yet an another temperature, but 
again by adding too many temperatures we lose our ability to ascertain the quality of the fit. 
Notably, A^(8) is also underestimated by the fit for HD 38087, and A^(7) is underestimated 
by the fits for HD 149404, HD 199579, and HD 210839 (the underestimation in smallest 
in the last case) — all cases of the column density of the highest detected J level being 
underestimated. A second possibility is that our column densities for the limiting cases are 
overestimated due to random noise enhancements of the relatively few detected lines (though 
we note we are confident that the detections are real due to consistency with velocity offsets). 

The weighted average value of T02 is 77 K, which is higher than the average Tqi de- 
termined in the FTCS — reasons for this discrepancy were discuss ed above. The weighted 



average of Ts4- is 2 78 K, which is consistent with previous surveys (jSpitzer fc Cochran 



1973 



Spitzer et al.lll974h . 



4.3. Meudon PDR Code Modeling 



As no ted in TO iNehme et al.l (|2008[ ) used the Meudon PDR code (most recently de- 



scribed by Le Petit et al. 20061 ) to model the HD 102065 line of sight. They were able to 



replicate certain observational properties of this line of sight, but not the higher- J H2 or 
CH"*" column densities. Increasing the radiation improved the fit to these excited column 
densities, but resulted in a worse overall fit to all of the observational quantities. 

We also used the Meudon PDR code to model two of our lines of sight. We take 
an approach similar to iNehme et al.l (120081 ): we fix all the possible inputs {Ay, Ry, total 
hydrogen column density, etc.) and allow only the density and radiation field to vary. Our 
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grid is very coarse, but covers a wide range in each variable. The local density nn is allowed to 
vary from 20 — 10^ cm~^, and values of the radiation field are allowed to vary from 0.14 — 140 
G (where G is the average Galactic value assumed by the code). The spacing is logarithmic, 
by factors of ~ 1.5 (~ 0.15 dex on the logarithmic scale). 

We selected two lines of sight, HD 46056 and HD 147888, as test cases. Neither line of 
sight necessarily represents a "typical" line of sight in our target list. It would be incorrect 
to say that such a line of sight exists in our sample, as many of the lines of sight were chosen 
for the FTCS due to their interesting properties. HD 46056 is in the open cluster NGC 2244 
and has a slightly below-average value of Ry and the subsequently expected higher far-UV 
extinction. In addition, the 2175 A extinction feature is weaker and narrower than in other 
lines of sight, qualities often associated with steep far-UV extinction. HD 147888 is notable 
as a sight line that passes through the p Oph complex]^ and at ~ 130 pc has one of the 
shorter line of sight pathlengths in the FTCS. HD 147888 has a very large value of Ry (> 4), 
and therefore less far-UV extinction. The fixed parameters for these lines of sight are listed 
in Table [H 

We perform the calculation using logarithmic column densities and errors of both 
H2 and the other molecules. This is done for two reasons. First, it weights the measured 
column densities more evenly, in better proportion to the true uncertainties. Second, when 
the calculation is performed using linear column densities, the best fit corresponds to a point 
at the density edge of our grid (and near the radiation edge) in both lines of sight. This was 
true for an earlier smaller and even more coarse grid that we calculated, but expanding the 
grid and increasing its resolution did not resolve the issue. However, a local minimum in the 
contour plots was found elsewhere in the grid for both lines of sight. These local minima 
are similar to the minima obtained when considering logarithmic column densities. 

The solutions for both lines of sight, reported in Table [12] are very high-density, h igh- 



radiation solutions. This is quite dissimilar to the results found by iNehme et al.l ( 120081 ) for 
HD 102065. We note that our solutions to both hues of sight are poor overall, with ^ 1, 
although some of the H2 column densities are nearly within the errors for HD 147888. In 
particular, the J = — 1 column densities, molecular fraction of hydrogen, and the ortho- 
to-parahydrogen ratio are very poorly fit. Our best f it mod els also severely underpredict the 



CH"*" column densities, similar to the iNehme et al.l (120081 ) models, and overpredict the CH 
column densities. 



Furthermore, with regard to HD 147888, we note that ISonnentrucker et al.l (j2007| ) also 



made an estimation of the density for this line of sight from C2 excitation that is significantly 



i^HD 147888 is also known as p Oph D. 
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smaller (215 cm~^) than our solution. This result assumed a typical interstellar radiation 
field, although lower K I and Na I abundances in this line of sight also imply an above 
average radiation field which our results do imply. 



Both our results and the iNehme et al.l ( l2008l ) results show the d ifiiculty of applying the 
Meudon PDR code to the issue of H2 excitation. The quality of the iNehme et al.l ( l2008l ) fit 
is much better than ours — their result for HD 102065 had a value of xl that was less than 
1, while our results for the xl of the best fit are significantly greater than 1. 

One important shortcoming of the Meudon PDR code, which is also true of many 
radiative transfer codes, is that such codes are usually steady-state models. However, the 
ISM is quite dynamic and r arely in comp l ete lo cal thermodynamic equilibrium. Specifically, 
shocks are not considered. iNehme et al.l ( l2008l ) concluded that shocks or other dissipation 
of turbulent thermal energy are the best explanation for the elevated levels of CH"*" and 
higher-J H2 relative to their best fit model's predictions. 

This same concern applies to our results as well. However, it is only one of several possi- 
ble explanations for the discrepancy of our measurements with the results of our test models. 
These include the following, listed roughly from greater to lesser probable importance: 



1. The steady-state assumption of the Meudon PDR code, especially the lack of shocks 

2. The assumption of a single, uniformly- dense cloud as opposed to multiple cloud com- 
ponents of varying density 

3. Errors in the input measurements, e.g. the column density of H I (particularly for HD 
147888, which has an unusually high value of iV(H I) compared to the other line of 
sight parameters and is uncertain due to the late spectral type of the star) 

4. Possible systematic errors have not been quantitatively accounted for, which would at 
least reduce the 

5. Chemistry and dust parameters were not tweaked, but were left at default values — 
although the values are reasonable and consistent with Galactic averages 

6. The coarseness of our grid makes any solution less than optimal 



We note that a different radiative transfer code by iBrowning et al.l (120031 ) was able 
to produce much better results for the higher-J H2 using either high radiation fields or 
"concatenated" models, that is, assuming a s uperposition of c l ouds as seen by the observer 



However, one major difference between the iBrowning et al.l (120031 ) and the iNehme et al. 
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( 120081 ) is that the former requires the input of both a density and a temperature which 
are assumed to be uniform throughout the cloud, while the latter more realistically only 
assumes a uniform density (which still may be a poor repre sentation) and calculates the 
temperature throughout the cloud. The iBrowning et al.l ( 120031 ) fits also did not consider the 
other observational data that is available that, in principle, any adequate model should be 
able to match. 



It remains to be seen whether a concatenation in the manner of iBrowning et al.l (j2003[ ) 
can be used in conjunction with a more complex cod e to consistently produce better results 
for a wide range of lines of sight. iNehme et al.l ( 120081 ) do attempt this for their line of sight, 
but find that it requires a second cloud of warm (~ 250 K), dense (uh > 10^ cm~^) gas within 
0.03 p c from the ba c kgrou nd star that has a ortho-to-para hydrogen ratio at formation near 
unity. iNehme et al.l (|2008[ ) note, in addition to the extreme nature of these conditions, that 
IRAS observations imply a distance of the cloud at least 0.12 pc from the star and an ortho 
to-par a formation ratio of ~ 1 is theortetically and experimentally unfounded. iNehme et al. 



(120081 ) therefore conclude that this scenario is not plausible, and favor the explanation of 
the dissipation of thermal energy to explain the higher- J and CH"^ results. We further note 
that given the qualitative similarity of our column density results for other lines of sight, 
this solution is even less plausible as a widespread explanation for the observations. 

We note that through this work we have substanti ally increased the h igher- J H2 mea- 
surements for lines of sight with N{B.2) > 10^° (see iNehme et al.l |2008| . and references 
therein). This provides a potential basis for additional modeling work that could explore 
questions such as whether a concatenation model is potentially a feasible explanation for 
certain lines of sight. 



5. SUMMARY 

We have undertaken a study of 22 lines of sight with properties that qualify them as 
"translucent lines of sight" though not necessarily "translucent clouds" , per the distinctions 



discussed in the FTCS and lSnow fc McCalll ( l2006l ). We find, consistent with previous studies 



of the interstellar medium, that H2 in the J > 2 states is overabundant relative to the 
expectation of a single, thermalized cloud. We do not find conclusive evidence that the 
observed lower- J and higher- J H2 are dominated by physically distinct clouds in any of 
these lines of sight — which we searched for in the form of trends in velocity offset and/or 
6-value as a function of J. We are unable to find a convincing correlation, either positive 
or negative, between 6-value and J-level. Similarly, we do not see significant evidence that 
there are trends between velocity offsets of the H2 transitions and J. However, the results 
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for both 6-values and velocity offsets are subject to tlie systematic errors discussed in §3.2[ 
and we are unable to state a null result more strongly. 

We explored the possibility that the column densities and abundances of various atomic 
and molecular species are correlated with the column densities of the different H2 J-levels 
and interesting derived ratios such as A^(4)/A^(0). We find many correlations with a high 
calculated confidence level, and virtually all of these have a straightforward physical inter- 
pretation. These correlations exist even when accounting for the common correlation with 
total column density. 

Further modeling is needed to understand the J > 2 H2 in lines of sight such as these. 
Our limited modeling attempts using the Meudon PDR code with two lines of sight resulted 
in solutions that are a poor match to several observed quantities such as the column densities 
of CH, CH+, ^^CO. The deviation of these lines of sight from their representation as single 
clouds in a steady state (particularly the influence of shocks) are likely to be the main reasons 
for the lack of quality solutions. 
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fits in addition to the other atomic and molecular column densities that are available on his 
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HD 38087 
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Fig. 1. — Above is an interesting region of the FUSE spectrum of HD 38087, showing H2 
hnes up to J = 8. The tick marks are at laboratory wavelengths; in this spectrum, the lines 
are shifted by ~ 11 km s^^ (< 0.5 A). J levels are in parentheses. HD 38087 shows the only 
clear detections of J = 8 in our sample of 22 lines of sight; the detected J = 8 line here is 
bluer feature of the two features near 1042.8 A. 
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Fig. 2. — Column density results of several different variations of the curve-of-growth method 
for HD 38087. Brief color-coded descriptions are given in the plot; see text for further 
details. Slight offsets along the abscissa are introduced for clarity. iV(0) and A^(l) are not 
independently derived, and thus are the same for all methods. 



- 36 - 



Curves of Growth, HD 24534 through HD 185418 
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Fig. 3. — Curve of growth fits for the HD 24534 through HD 185418 hues of sight, assuming 
a uniform 6-value for all J levels. The error bars of the data points are color-coded by J 
level. The black lines are single component curves of growth with the 6-value of the best 
fit solution. The divergence of the lines to the right represent the range in values of the 
damping constant 7, which is different for each transition. 3-0" limits of the first undetected 
J- level and the J = 0, = 1 vibrational level are represented with black triangles and 
purple squares (with downward lines), respectively; detected vibrational lines are blue-green 
circles (J = 0) and green diamonds (J = 1), with matching colored error bars. Purple: — 
J = 2. Blue: — J = 3. Light blue-green: — J = 4. Green: — J = 5. Yellow: — J = 6. 
Orange: — J = 7. Red: — J = 8. 
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Curves of Growth, HD 186994 through HD 210839 




log(NfA) (NfA in cm"') 



Fig. 4.— Same as Fig. [3l except for the HD 186994 through HD 210839 hnes of sight. 
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Curve of Growth with Assumed Velocity Structure, HD 210839 
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Fig. 5. — Same as Fig. [3l except the specific case of HD 210839 when fixed to the velocity 
structure of CH. The equivalent widths fit to the curve of growth in this line of sight much 
better with a much higher 6- value than the velocity structure assumed here. The results are 
similar whether CH or CH"*" data is used, and similar effects are observed in some other lines 
of sight, particularly HD 199579. 
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Fig. 6. — Left: — The distribution of the calculated rank correlation coefficients of 6-values 
vs. Ej/k. The abscissa is the significance of the rank correlation coefficients converted into 
a Gaussian cutoff value. Right: — The distribution of slope significance (slope divided by 
slope error) for the same data set. Note that positive Gaussian cutoff values correspond to 
negative correlations/slopes. Also note that there are 68 data points — 17 lines of sight with 
four data segments each. Reasons that there are apparent anticorrelations are discussed in 
the text. 
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Fig. 7. — Same as Fig. [6l except using 6-values from curve-of-growth method. For this 
reason, there are only 22 data points. 
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Fig. 8. — Same as Fig. [6l except for velocity offsets. 
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Fig. 9. — Fits of two excitation temperatures to statistically weighted column densities 
for the HD 24534 through HD 185418 lines of sight. Errors are exaggerated on the plot's 
logarithmic scale by a factor of 5 for clarity. Filled circles are the main J = — 8, z/ = 
states. The 3-cr upper limit on the first undetected J-level is a downward arrow with a 
horizontal bar. 3-0" upper limits on the J = 0, u = 1 state are shown with downward arrows 
without a bar. Detections of vibrationally excited material are shown with open squares, 
only applicable for HD 38087 on this plot and HD 199579 in Fig. [TOl See text for further 
discussion. 
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Fig. 10.— Same as Fig. |9] except for the HD 186994 through HD 210839 hnes of sight. 
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Table 1. Lines of Sight: Stellar Data 



Star Name 


Spectral Class 


1 


b 


Distance (pc) 


Distance Ref. 


HD 24534 


09.5pe 


163.08 


-17.14 


590 


1 


HD 27778 


B3V 


172.76 


-17.39 


223 


2 


HD 38087 


B5V 


207.07 


-16.26 


480 


1 


HD 40893 


BOIV 


180.09 


+4.34 


2800 


1 


HD 46056 


08V 


206.34 


-2.25 


2300 


1 


HD 46202 


09V 


206.31 


-2.00 


2000 


1 


HD 53367 


BOIV:e 


223.71 


-1.90 


780 


1 


HD 147888 


B3/B4V 


353.65 


-fl7.71 


136 


2 


HD 149404 


09Ia 


340.54 


+3.01 


820 


1 


HD 170740 


B2V 


21.06 


-0.53 


213 


2 


HD 179406 


B3V 


28.23 


-8.31 


160 


1 


HD 185418 


B0.5V 


53.60 


-2.17 


950 


1 


HD 186994 


BOHI 


78.62 


+10.06 


2500 


1 


HD 192639 


07.5IHF 


74.90 


+1.48 


1100 


1 


HD 195965 


BOV 


85.71 


+5.00 


1300 


1 


HD 197512 


BIV 


87.89 


+4.63 






HD 199579 


Oe.SHI 


85.70 


-0.30 


1100 


3 


HD 203938 


B0.5IV 


90.56 


-2.23 


700 


1 


HD 206267 


06(F) 


99.29 


+3.74 


1000 


4 


HD 207198 


09n 


103.14 


+6.99 


1000 


1 


HD 207538 


09.5V 


101.60 


+4.67 


880 


1 


HD 210839 


06e 


103.83 


+2.61 


505 


2 



References. — (1) From the DIBS database, |http:/ /d ib. uiuc.edu/', values com- 
piled by L. M. Hobbs, spectr oscopic distance mo dulus. (2) Hipparcos value of 4- 
a significance or better. (3) ISavaee et al.l lll985l ). (4) From the DIBS database, 
http://dib.uiuc.edu/, values compiled by L. M. Hobbs, stellar association. 
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Table 2. Lines of Sight: Reddening Data 



Sightline 


Eb^v 


Ref. 


Av 


Ref. 


Rv 


Ref. 


HD 24534 


0.59 


1 


2. 


.05 


2 


3.47 


2 


HD 27778 


0.37 


3 


1. 


.01 


2 


2.72 


2 


HD 38087 


0.29 


3 


1. 


.61 


4 


5.57 


4 


HD 40893 


0.46 


3 


1. 


.13 


4 


2.46 


4 


HD 46056 


0.50 


3 


1. 


.30 


4 


2.60 


4 


HD 46202 


0.49 


3 


1. 


.39 


4 


2.83 


4 


HD 53367 


0.74 


3 


1. 


.76 


4 


2.38 


4 


HD 147888 


0.47 


3 


1. 


.91 


4 


4.06 


4 


HD 149404 


0.62 


5 


2. 


.19 


5 


3.53 


5 


HD 170740 


0.48 


3 


1. 


.30 


2 


2.71 


2 


HD 179406 


0.33 


3 


0, 


.94 


4 


2.86 


4 


HD 185418 


0.50 


3 


1, 


.16 


2 


2.32 


2 


HD 186994 


0.17 


3 


0, 


.53 


6 


3.10 


7 


HD 192639 


0.66 


1 


1, 


.87 


2 


2.84 


2 


HD 195965 


0.25 


1 


0, 


.77 


8 


3.08 


6 


HD 197512 


0.32 


9 


0, 


.75 


2 


2.35 


2 


HD 199579 


0.37 


9 


1, 


.09 


2 


2.95 


2 


HD 203938 


0.74 


3 


2, 


.15 


2 


2.91 


2 


HD 206267 


0.53 


3 


1, 


.41 


2 


2.67 


2 


HD 207198 


0.62 


3 


1, 


.50 


2 


2.42 


2 


HD 207538 


0.64 


3 


1, 


.44 


2 


2.25 


2 


HD 210839 


0.57 


1 


1, 


.58 


2 


2.78 


2 



References. — (1) iDiplas fc Savagj l ll994l) . (2) 
iRachford et al.l ll2002bl ): if Ay, reference 6 also applies. (3) 
From the DIBS database, http://dib.u iuc.edu/ , values by the 
DIBS team. (4) IRachford et al.l ||2009| '): if Ay, reference 6 also 
applies. (5) IValencic et al. I 1I2OO4I '). (6) Calculated through 
the definition Ry = Ay/Eg_y using values from the other 
two columns. (7) Ry not independently measu red and is as- 
sumed to be the Galactic average of 3.1. (8) iNeckel et al] 
il98di) . (9) lAiello et"an lll988l'l . 



Table 3. Measured Equivalent Widths 



Wavelength 


J Level 


HD 24534 


HD 27778 


HD 38087 


HD 40893 


HD 46056 


HD 46202 


HD 53367 


HD 147888 


HD 149404 


HD 170740 


Line of Sight 
HD 179406 HD 185418 


t 


919.048 


4 




























920.241 


2 








118.8 it 53.5 




128.2 it 28.1 






166.4 it 58.0 






104.4 it 29.6 


i 


921.730 


3 
























87.3 it 25.5 




922.303 


4 




























922.893 


3 




























927.017 


2 


46.1 it 12.2 






166. ± 19.0 




161.3 zb 43.9 






136.5 zb 17.3 






107.5 zb 26.0 




928.436 


3 








142.1 zb 23.6 


114.3 it 44.3 


97.3 zb 36.9 






128.7 it 33.7 






51.4 ± 21.4 




932.604 


2 








160.8 it 27.9 




123.3 it 28.6 






147.5 it 39.5 






101.6 it 20.3 




933.578 


3 








175.1 it 17.3 


106.9 it 35.1 


183.9 it 54.1 






120.9 it 16.3 






93.3 it 12.7 


] 


933.788 


4 










82.8 zb 31.7 


92.0 zb 33.7 






94.6 it 12.1 






40.6 ± 12.3 




934.141 


2 


61.0 ± 8.9 






125.6 zb 17.7 


99.6 zb 21.8 


151.0 zb 43.8 


175.0 zb 72.9 




109.5 ± 14.5 






93.6 zb 9.5 




934.789 


3 


47.9 ± 9.8 






158.4 ± 35.6 


108.8 ± 26.8 


131.3 ± 19.1 


67.5 ± 35.8 




122.1 zb 25.7 






71.8 ± 12.7 


i 


935.958 


4 


25.8 it 7.8 






80.0 it 12.5 




72.1 it 15.8 






93.7 zb 35.7 






50.1 it 14.5 




936.463 


5 


















68.3 zb 9.1 










936.854 


3 




























937.551 


4 




























938.726 


4 




























938.909 


5 




























940.384 


4 


















43.0 it 12.0 






17.2 ± 15.4 




940.623 


2 








161.5 it 24.6 


123.3 it 22.8 


168.2 it 37.4 






141.6 it 12.5 


226.1 it 106.4 


\''\ 


112.5 ± 8.3 


I 


940.882 


5 












94.4 ± 25.7 






59.1 it 13.1 




. r\ 

■ ■ 




2 


941.596 


2 








171.0 ± 41.3 


111.0 i 24.6 


151.3 ± 38.3 


222.0 it 74.5 


69.6 ± 33.4 


142.7 ± 9.8 


191.9 it 84.9 


. . 


105.3 ± 10.0 




942.685 


5 










36.0 it 34.2 


52.8 it 27.8 






51.9 dt 12.4 


64.8 ± 56.1 


... 1 






942.962 


3 


35.4 ± 10.0 


96.6 ± 62.7 


144.0 it 39.7 


153.6 it 18.3 


94.6 it 22.8 


109.0 zb 20.8 




95.1 ± 32.8 


109.2 it 11.8 






78.4 it 10.6 




944.328 


3 


32.9 ± 10.1 






145.5 it 17.1 


126.9 it 28.4 


119.6 zb 22.1 




75.8 ± 26.4 


127.2 it 17.0 


54.9 dt 21.0 




65.1 it 11.2 


( 


944.718 


5 












47.7 ± 29.6 
















947.885 


4 




























951.672 


3 


76.8 ± 8.7 






173.8 ± 24.0 


127.7 it 20.8 


138.0 it 18.0 


51.7 ± 28.2 




134.7 it 13.5 






88.0 ±8.1 


■ 


955.851 


4 


























2 


956.578 


2 


124.7 it 7.5 


169.2 it 46.5 


224.5 dt 46.4 


182.0 it 20.7 


137.3 it 21.7 


172.3 it 20.6 




147.9 it 18.0 


187.4 it 34.9 


163.3 it 36.9 




137.1 it 22.1 




957.406 


6 






14.6 it 10.5 






















957.650 


2 


102.7 ± 5.2 


148.8 ± 21.9 






130.8 zb 24.4 


174.1 zb 16.5 




132.3 zb 14.0 


145.6 zb 31.5 




126.5 ± 53.9 


105.6 zb 8.2 




957.818 


5 










96.7 it 44.0 


64.6 it 15.6 




35.9 it 12.6 


80.6 it 8.4 






20.7 ± 8.8 




958.009 


5 










60.4 ± 35.9 


43.3 ± 15.7 




29.3 ± 13.8 


65.6 ± 13.3 


31.2 it 18.0 




32.0 ± 13.6 




958.945 


3 


53.7 ± 9.3 


92.6 dt 21.1 


102.5 it 16.7 


152.5 ± 22.6 


149.4 ± 25.9 


104.9 ± 16.4 


77.8 ± 42.2 


101.9 ± 20.3 


129.3 ± 15.0 






88.5 ± 10.0 


i 


960.263 


5 












40.9 ± 20.8 






49.4 it 12.7 










960.449 


3 


85.2 it 7.5 




73.9 it 18.2 


141.3 it 20.8 


126.6 it 20.0 


141.7 it 14.6 


107.6 it 54.1 


55.3 it 19.8 


113.7 it 15.7 










962.151 


4 






42.9 it 23.9 


71.8 ± 16.2 


95.7 ± 42.9 


80.7 it 16.9 




24.6 ± 21.2 


117.2 zb 16.6 


86.9 ± 31.0 




62.9 zb 11.3 




965.791 


2 


























1 


966.779 


3 


















137.5 ± 28.3 






82.6 lb 32.2 


i 


967.278 


2 


207.8 it 15.6 






221.1 ± 35.1 


186.0 it 51.6 


234.8 it 70.6 




287.3 ± 43.2 


226.3 ± 31.1 






221.4 ± 48.1 


1 


967.673 


3 


53.4 it 8.9 






157.9 it 31.7 




212.0 it 41.5 




69.5 ± 19.7 


102.0 it 15.7 






87.0 it 10.1 


■ 


968.292 


2 


145.8 it 9.1 






175.9 it 13.0 


144.3 it 29.5 


180.0 it 30.3 


375.6 it 73.3 


185.8 it 26.0 


170.2 it 68.3 


283.0 it 80.5 




129.0 ± 9.1 


1 


968.555 


5 




























968.664 


4 
























37.6 ± 8.1 




970.560 


3 


37.6 ± 10.3 






139.7 ± :n.7 


73.8 zb 21.3 


138.8 zb 47.6 












68.7 ± 8.3 


i 


970.835 


4 


29.4 ± 8.1 








53.0 ± 17.2 


71.4 ± 21.2 












45.1 rt 10.5 




971.074 


5 








43.9 it 19.4 
















21.0 it 11.0 




971.387 


4 




























974.156 


2 


207.8 it 9.6 






















217.0 it 25.1 


1 



Table 3 — Continued 



AVavelengtn 


J Level 


ilU J40o4 






nU 4Uoyo 


iiu 4dUod 


iiu 4oJUJ 


nL) OOOD / 


iiu i4/ocSo 


XJT^ ^ ACtACiA 

ilU i4y4U4 


nL) 1 / U / 4U 


Line 
iiu i / y4uo 


of Sight 

iiu loo41o 




974.286 


5 
























50.8 ± 23.9 




974.884 












47.9 ± 18.0 


71.7 ± 12.8 






74.7 lb 10.0 






28.4 ± 12.9 




97.S.344 


2 


107.2 ± 5.8 






207.4 ± 36.4 


90.5 ± 19.6 


157.5 ± 22.0 


425.8 ± 180.1 




155.7 ± 22.2 


218.2 ± 76.3 




130.2 ± 13.9 




978.217 


3 


56.0 ± 6.4 




92.2 ± 35.8 


170.3 ± 26.2 


137.6 ± 25.6 


176.4 ± 19.1 




48.1 ± 27.4 


118.5 ± 17.1 






76.9 ± 15.7 




979.803 


4 








81.4 ± 18.6 


84.4 ± 22.3 








113.1 ± 7.5 










980.498 


6 






62.4 ± 26.0 






















983.589 


2 


162.3 ± 11.8 






216.7 ± 27.3 


135.3 ± 33.0 


163.2 ± 22.9 




253.0 ± 20.2 








185.1 ± 23.4 




984.862 


2 




























986.241 


2 


























1 


987.44."} 


3 




























987.767 


3 
























115.2 ± 21.7 




987.972 


2 
























272.9 ± 33.2 




989.086 


2 












248.1 ± 15.0 




135.2 ± 16.6 


145.1 ± 25.1 






134.6 ± 11.9 




989.553 


4 
























54.7 ± 16.4 




994.227 


4 


22.6 ± 6.6 


63.6 ± 23.9 




77.4 ± 14.3 


98.3 ± 16.6 


122.1 ± 17.9 


112.8 ± 38.6 


64.0 ± 14.4 


119.6 ± 22.2 


88.0 ± 26.0 


71.7 ± 20.3 


47.4 ± 8.5 




994.517 


7 






22.4 ± 20.3 






















994.871 


2 


179.8 ± 7.7 


























994.924 


5 




























997.638 


5 




46.5 ± 21.1 


51.6 ± 23.8 




67.4 ± 19.0 


58.2 ± 14.6 


54.6 ± 30.1 




86.7 ± 16.6 






10.3 ± 9.6 




997.824 


3 


52.3 ± 4.6 


84.4 lb 12.9 


98.8 ± 20.2 




120.3 lb 13.9 


142.1 ± 16.0 


104.1 ± 26.4 


61.4 ± 9.7 


122.8 ± 17.3 


74.7 ± 13.1 


78.8 ± l3.8 


80.1 lb 6.8 




998.331 


6 






49.1 ± 10.6 
















■ ■ ir\ 






998.425 


6 






27.1 ± 7.2 












22.8 ± 8.3 




. . --J 






999.268 


4 


:-»;-!.() ± 4.4 


40.8 ± 11.6 


80.0 ± 10.5 


67.8 ± 7.8 


80.5 ± 12.2 


80.4 ± 9.3 


55.6 ± 20.1 


65.1 ± 9.9 


109.4 ± 14.7 


74.4 ± 20.5 


57.0 ± p.l 


54.8 ± 5.6 




1003.982 


2 


207.5 ± 7.4 


216.1 ± 21.8 


438.0 ± 66.3 


230.0 ± 11.9 


218.9 ± 18.6 


269.0 ± 24.6 




262.4 ± 28.3 


228.6 ± 25.6 


294.9 ± 31.5 




181.5 ± 7.4 




1005.390 


2 


157.1 lb 5.3 


235.6 ± 16.8 


456.9 lb 32.4 


219.1 ± 13.9 


183.6 lb 14.3 


200.1 ± 42.7 


433.7 ± 71.7 


291.2 ± 11.9 


222.4 ± 21.3 


310.1 lb 28.7 


161.6 ± 10.6 


161.0 ± 5.6 




1005.802 


7 






33.3 ± 9.5 






















1009.024 


2 




























1010.129 


3 




























1010.938 


2 


243.6 ± 7.7 


227.6 ± 39.4 


509.7 ± 86.6 


238.1 ± 12.3 


245.9 ± 26.2 


209.5 ± 30.9 




301.7 ± 26.6 


275.8 ± 41.6 


379.6 ± 36.9 


177.3 ± 13.6 


224.6 ± 8.7 




1012.169 


2 




























1012.259 


4 




























1014.504 


3 




























1016.458 


2 


184.6 ± 3.9 


182.3 ± 8.3 


349.4 ± 22.7 


200.9 ± 11.8 


174.4 ± 9.4 


229.2 ± 13.0 


392.8 ± 23.2 


224.0 ± 9.9 


205.6 ± 16.0 


316.6 ± 26.8 


164.9 ± 16.6 


166.6 ± 4.7 




1016.742 


6 






21.9 ± 8.6 






















1017.001 


5 


5.9 ±3.1 




51.1 ± 10.9 


31.6 ± 7.3 


66.8 ± 10.1 


48.8 ± 7.7 


46.2 ± 28.6 


49.6 ± 13.0 


75.8 ± 7.2 


36.9 lb 9.5 




20.0 lb 4.4 




1017.831 


5 


22.3 ± 3.8 






48.6 ± 5.9 


67.8 ± 9.6 


96.3 ± 6.9 




52.0 ± 7.2 


106.6 ± 3.9 


38.0 ± 12.9 


53.4 ± 20.8 


30.6 ± 3.9 




1018.214 


7 






19.6 ± 11.5 






















1019.013 


6 






37.0 ± 9.5 








17.3 ± 8.5 




22.6 ± 4.0 


16.9 ± 9.9 








1019.500 


3 


52.4 ± 2.9 


89.0 ± 12.7 


154.7 ± 13.3 


187.2 ± 15.9 


112.4 ± 8.7 


166.2 ± 10.4 


117.7 ± 10.4 


69.2 ± 8.6 


134.9 ± 15.7 


97.7 ± 11.3 


136.6 ± 27.0 


90.1 ±4.0 




1021.207 


6 






47.3 ± 13.7 






















1022.580 


6 


















33.6 ± 8.1 










1023.434 


4 


38.8 ± 8.9 






89.0 ± 16.5 


82.2 ± 21.4 


117.0 ± 19.4 












38.6 ± 8.2 




1026.526 


2 


























1 


1028.801 


7 




























1028.985 


3 


75.7 ± 3.1 


89.1 ± 21.0 


192.3 ± 26.6 


187.4 ± 9.8 


135.8 ± 14.2 


183.0 ± 20.1 


109.6 ± 9.9 


84.6 ± 25.6 


152.6 ± 24.0 


88.0 ± 12.2 


94.6 ± 13.9 


91.9 ± 4.4 




1030.071 


6 






57.7 ± 18.3 




17.7 ± 16.1 








28.9 ± 13.6 










1031.191 


3 


69.1 ± 2.7 


81.3 ± 13.5 


178.3 ± 22.5 


172.2 ± 13.1 


120.9 ± 8.6 


153.9 ± 8.8 


104.6 ± 9.9 


89.9 ± 10.4 


135.2 ± 18.6 


90.0 ± 11.7 


112.6 ± 9.2 


89.2 ± 4.0 




1032.349 


4 


37.0 ± 2.8 


62.6 ± 8.9 


86.2 ± 11.5 


75.4 ± 5.9 


81.9 ± 8.4 


94.9 ± 5.9 


81.2 ± 8.2 


46.9 ± 7.2 


110.2 ± 12.3 


80.5 ± 8.8 


122.8 ± 63.9 


53.1 ± 3.1 




1033.914 


6 






60.8 ± 14.6 












17.5 ± 5.2 


51.1 ± 24.1 








1035.181 


4 


26.7 ± 4.2 


49.1 ± 13.8 


69.6 ± 12.7 


84.3 ± 7.4 


86.1 ± 11.1 


104.3 ± 11.4 


107.8 ± 22.8 


44.6 ± 7.7 


107.9 ± 9.4 


103.6 ± 14.2 


56.1 ± 12.9 


52.4 ± 4.8 





Table 3 — Continued 



Wavelength 


J Level 




nu J / / /8 


nL) ooUo / 


nL) 4Uoyo 


iiu 4dUod 






iiu 14(ooo 




nU 1 (U / 4U 


Line of Sight 
HD 179406 HD 185418 


HI 


1038 


689 


2 




276.1 ± 20.0 


423.9 ± 77.8 




234.5 ± 15.8 


209.3 ± 26.0 




273.0 ± 6.9 






189.2 ± 16.4 


246.4 ± 9.3 


12 


1040 


058 




15.3 ± 2.9 




77.8 ± 9.6 


17.3 ± 4.5 


43.7 ± 7.5 


72.7 ± 7.8 




28.2 ± 5.0 


93.5 ± 4.5 


53.1 ± 9.1 


35.3 ± 22.6 


46.1 ± 5.0 


i: 


1040 


366 


2 


200.1 ± 3.0 


209.5 ± 7.5 


354.7 ± 17.4 


212.3 ± 10.0 


189.2 ± 7.5 


221.1 ± 15.4 


410.7 ± 72.3 


217.3 ± 6.4 


267.4 ± 18.8 


319.7 ± 16.8 


246.0 ± 27.1 


197.1 ±4.1 


11 


1041 


158 


3 


63.0 ± 2.5 


103.9 ± 6.4 


219.4 ±8.1 


182.1 ± 7.2 


128.6 ± 6.9 


197.9 ± 17.8 


141.8 ±7.7 


82.8 ± 6.6 


140.2 ± 8.4 


125.2 ± 7.6 


107.6 ± 16.0 


99.3 ± 3.6 


10 


1041 


729 


6 






42.8 ± 6.4 








15.6 ± 8.0 




15.9 ± 3.6 


35.2 ±8.1 


5.4 ± 4.2 






1042 


742 


8 






16.2 ± 6.0 






















1043 


502 


3 


61.2 ± 2.6 


89.7 ± 4.5 


170.9 ± 6.9 


178.7 ± 8.8 


123.5 ± 7.1 


159.2 ± 11.8 


111.1 ± 6.4 


64.2 ± 4.6 


133.2 ± 6.4 


110.4 ± 6.1 


129.1 ± 16.8 


89.8 ± 2.9 


9' 


1044 


542 


4 


36.6 ± 2.6 


49.0 ± 5.3 


75.1 ± 6.1 


102.1 ± 11.8 


79.4 ± 7.6 


88.6 ± 4.8 


74.6 ± 6.5 


52.2 ± 4.8 


109.7 ± 10.6 


60.0 ± 5.6 


81.6 ± 21.1 


52.9 ± 3.3 


5( 


1045 


801 


6 






37.1 ± 5.6 










15.7 ± 5.3 


11.9 ± 3.0 










1047 


550 


4 


32.3 ± 3.3 


41.9 ± 5.7 


63.7 ± 6.9 


79.3 ± 7.4 


77.9 ± 8.7 


85.5 ± 6.5 


68.0 ± 8.8 


36.2 ± 5.0 


105.5 ± 8.6 


61.6 ± 8.0 


51.9 ± 4.8 


45.9 ± 3.4 


4^ 


1047 


711 


7 






28.7 ± 7.4 






















1048 


829 




























2 


1051 


498 


2 




338.0 ± 19.8 


406.3 ± 100.1 




288.4 ± 29.7 






262.4 ± 6.7 






137.7 ± 9.2 


207.7 ± 17.6 


12 


1052 


393 


7 






21.6 ± 4.7 






















1052 


496 


5 


9.0 ± 2.9 


9.0 ± 7.2 


47.8 ± 8.4 


27.4 ± 6.1 


62.2 ± 8.0 


61.7 ± 6.7 


45.6 ± 7.4 


29.4 ± 4.7 


91.9 ± 4.1 


46.4 ± 8.3 


20.2 ± 6.0 


22.5 ± 4.0 


1' 


1053 


976 


3 


66.7 it 2.3 


91.1 it 4.7 


190.3 ± 6.9 


179.3 ± 8.7 


129.2 ± 6.6 


173.6 ± 12.4 


125.6 ± 6.7 


78.9 ± 4.7 


150.6 ± 10.1 


114.9 ± 6.8 


105.2 ± 16.4 


98.8 ± 2.7 


10 


1054 


515 


8 






4.1 ± 3.4 






















1056 


472 


3 


61.0 ± 2.5 


82.7 ±5.1 


143.1 ± 7.2 


179.4 ± 8.0 


133.7 ± 6.8 


161.9 ± 10.6 


124.5 ± 6.7 


62.6 ±4.1 


145.7 ± 16.6 


83.1 ± 4.9 


110.4 ± 12.8 


91.4 ± 2.7 


10 


1057 


380 


4 


33.8 ± 2.6 


40.3 ± 5.7 


72.6 ± 7.4 


89.6 ± 8.7 


86.2 ± 7.7 


92.9 ± 5.2 


73.5 ± 6.0 


47.8 ± 5.0 


113.9 ± 4.2 


66.5 ± 6.3 


60.2 ± 17.7 


47.3 ± 3.2 


5C 


1058 


316 


6 






36.6 ± 6.9 










12.2 ± 5.3 


13.6 ± 4.0 


14.5 ±8.1 


8.8 ± 3.|7 






1060 


031 


7 






32.1 it 7.3 












14.2 ± 4.4 










1060 


581 


4 


31.3 ± 3.0 


45.3 ± 6.3 


68.6 ± 7.6 


79.0 ± 9.2 


73.5 ± 7.3 


83.5 ± 7.6 


77.4 ± 10.0 


47.7 ± 5.2 


109.4 ± 11.4 


53.3 ± 7.1 


68.2 ± 


42.7 ± 3.4 


4-: 


lt)61 


697 




19.5 ± 4.4 


22.5 ± 9.0 


66.4 ± 9.6 


31.9 ± 6.7 


82.4 ± 8.9 


ta.'Z ± 7.4 


82.4 ± 18.8 


37.1 ± 7.5 


88.0 ± 8.3 




49.0 ± l!f.8 


26.8 ± 5.1 


21 


1064 


994 


2 




266.2 ± 44.7 






237.9 ± 33.3 


245.0 ± 33.7 




238.6 ± 9.9 






182.7 ±8.1 


248.0 ± 28.2 


11 


1065 


596 


5 






48.5 ± 9.2 


41.5 ± 6.1 


53.1 ± 9.1 


54.6 ± 6.6 


39.1 ± 6.4 


38.9 ± 4.7 


85.0 ± 15.0 


22.5 ± 6.1 


14.3 ± 4.7 


18.6 ± 4.7 


i; 


1066 


900 


2 


171.6 ± 2.6 


172.9 ± 6.9 


304.3 ± 16.6 




206.7 ± 8.6 


214.6 ± 14.6 




188.6 ± 6.3 


222.7 ± 26.9 


280.0 ± 9.4 


209.3 ± 19.6 


167.2 ± 3.7 


11 


1067 


479 


3 


60.6 ± 2.4 


99.0 ± 6.2 


172.2 ± 6.7 


179.1 ± 8.7 


139.9 ± 8.2 


163.4 ± 8.3 


116.6 ± 6.7 


71.1 ± 6.8 


166.7 ± 6.9 


99.2 ± 6.5 


104.4 ± 20.7 


87.9 ± 2.7 


9; 


1070 


141 


3 


60.6 ± 3.4 


96.3 ± 6.4 


161.7 ± 7.7 


184.2 ± 17.6 


127.6 ± 7.2 


162.1 ± 8.1 


104.6 ± 7.2 


69.6 ± 4.3 


136.9 ± 10.8 


95.2 ± 6.5 


127.7 ± 17.0 


91.0 ± 2.7 


9.. 


1070 


900 


4 


28.9 ± 3.5 


43.0 it 4.9 


59.8 ± 6.1 


118.6 ± 8.2 


87.7 ±8.1 


78.9 ± 6.9 


77.1 ± 9.1 


40.3 ± 4.6 


106.2 ± 6.9 


66.2 ± 6.6 


64.5 ± 13.2 


60.0 ± 3.2 


5^ 


1071 


497 


6 






27.4 ± 5.1 












14.9 ± 3.1 










1072 


501 


8 






15.1 ± 11.3 






















1072 


992 


7 






26.0 ± 6.0 






















1074 


312 


4 


26.7 ± 2.6 


63.3 it 6.2 


64.4 ± 6.2 


73.0 ± 7.6 


76.2 ± 8.1 


76.6 ± 6.4 


71.4 ± 7.1 


49.3 ± 6.8 


110.4 ± 4.1 


42.7 ± 6.5 


63.3 ± 17.0 


61.4 ± 3.4 


4S 


1075 


245 


5 


11.9 ± 3.4 


14.1 it 6.1 


61.3 ± 6.1 


42.7 ± 9.0 


46.7 ± 9.6 


54.6 ± 7.6 


37.6 ± 9.3 


26.2 ± 5.3 


81.1 ± 13.9 


39.3 ± 7.3 


32.2 ± 6.2 


23.7 ± 4.7 


IS 


1078 


266 


7 






29.2 ± 7.1 






















1079 


226 


2 
















174.7 ± 11.2 










10 


1079 


400 


5 










32.3 ± 11.0 


42.5 ± 16.1 




24.3 ± 5.9 








22.5 ± 8.9 


2' 


1079 


922 


8 






9.1 ± 4.4 






















1080 


492 


6 






28.6 ± 6.1 












13.7 ±3.4 


20.7 ± 7.6 


6.1 ± 3.3 






1081 


267 


2 


166.3 ± 3.9 


171.1 ± 6.6 


302.6 ± 9.7 


200.6 ± 17.2 


173.7 ± 9.6 


194.0 ± 7.6 


376.0 ± 9.8 


173.4 ± 6.2 


226.2 ± 12.3 


248.8 ± 8.7 


131.2 ± 17.7 


163.6 ± 4.6 


10 


1081 


713 


3 






166.0 ± 9.8 




99.2 ± 13.6 


161.3 ± 7.6 


112.4 ± 11.6 


60.7 ± 6.7 


134.0 ± 10.6 


74.6 ± 8.4 


120.6 ± 17.1 


87.5 ± 4.7 


9. 


1084 


562 


3 


92.5 ± 5.1 


76.7 ± 12.6 




214.6 ± 17.6 


169.6 ± 15.0 


200.9 ± 26.0 




91.4 ± 13.9 


174.7 ± 6.6 






84.5 ± 8.1 


11 


1085 


146 


4 


20.9 ± 5.9 


21.6 it 9.8 


66.8 ± 11.4 


64.5 ± 12.3 


67.3 ± 22.9 


66.9 ± 19.9 


76.8 ± 14.8 


31.9 ± 15.3 


110.3 ± 5.7 


33.2 ± 12.5 




25.6 ± 12.4 


42 


1085 


382 


6 






57.9 ± 21.8 






















1086 


630 


7 






35.2 ± 17.3 






















1088 


796 


4 


27.6 ± 3.6 


39.6 it 6.0 


69.1 ± 9.6 


49.6 ± 8.9 


60.8 ± 9.3 


64.6 ± 8.2 


64.7 ±8.1 


47.2 ± 6.2 


100.8 ±6.4 


46.9 ± 7.6 




36.1 ± 4.9 


3e 


1089 


515 


5 


7.8 ± 4.4 




59.6 ± 10.2 


29.1 ± 9.0 


44.4 ± 9.6 


54.4 ± 10.0 


36.0 ± 8.6 


31.6 ± 6.0 


67.0 ± 16.8 


23.7 ± 7.8 




15.4 ± 4.3 


6 


1094 


244 


2 
















139.7 ± 14.8 










10 


1096 


438 


2 


123.2 ± 4.3 


138.8 ± 6.3 


223.2 ± 8.4 




141.7 ± 7.0 


183.5 ± 6.3 


313.7 ± 6.4 


141.3 ± 5.2 


187.1 ± 16.3 


228.9 ± 7.6 


110.4 ± 20.1 


124.6 ± 3.5 


9' 


1096 


725 


3 


49.9 ± 5.1 


93.2 ± 9.1 


143.9 ± 13.9 




116.7 ± 10.6 


154.7 ± 7.0 


117.8 ± 8.8 


76.8 ± 8.7 


144.0 ± 13.5 




104.0 ± 11.1 


75.0 ± 4.1 


11 



Table 3 — Continued 



Wavelength 


J Level 


HD 24534 


HD 27778 


HD 38087 


HD 40893 


HD 46056 


HD 46202 


HD 53367 


HD 147888 


HD 149404 


HD 


170740 


Line of Sight 
HD 179406 HD 185418 


HD 186994 


HD 1 


1099.787 


3 


45.9 ± 3.3 


74.3 ± 6.2 


110.3 ± 7.3 


145.7 ± 7.1 


108.0 ± 7.2 


110.8 ± 6.2 


86.9 ± 8 


8 


60.7 ± 4.6 


128.5 ± 6.3 


82 


1 ± 7.6 


73.8 ± 9.1 


72.9 ± 3.3 


74.9 ± 4.5 


97.0 


1100.019 


6 






38.7 ± 8.2 




























1100.164 


4 


20.6 ± 3.3 


24.8 ± 6.1 


35.0 ± 6.5 


33.3 ± 8.0 


46.6 ± 11.2 


119.1 ± 7.8 


64.8 ± 6 


~y 


19.8 ± 3.6 


90.4 ± 10.9 


25 


± 5.5 


29.0 ± 9.3 


29.8 ± 3.8 


18.0 ± 3.6 


59.7 


1100.982 


7 






16.7 ± 3.7 




























1104.083 


4 




40.8 ± 2.9 


52.1 ± 3.2 


29.2 ± 3.7 


59.7 ± 4.6 


63.4 ± 3.9 


65.4 ± 3 


7 


34.5 ± 3.5 


92.0 ± 8.0 


61 


8 ± 5.0 


45.7 ± 8.5 


34.1 ± 2.6 


14.1 ± 4.0 


61.8 


1106.971 


7 






23.2 ± 4.4 




























1109.313 


5 






54.4 ± 4.7 




16.6 ± 4.6 


25.0 ± 4.1 


35.7 ± 4 


2 


31.0 ± 4.3 


47.1 ± 8.9 


27 


2 ± 5.2 




8.3 ± 2.2 




25.8 


1115.454 


6 






19.0 ± 4.1 




























1115.895 


3 


38.3 ± 1.5 


48.2 ± 2.3 


65.2 ± 4.0 




95.4 ± 4.9 


112.4 ± 3.4 


82.4 ± 3 


1 


42.8 ± 3.2 


134.1 ± 9.0 


53 


7 ± 4.7 


74.3 ± 6.2 


66.8 ± 1.9 


70.2 ± 2.2 


94.6 


1116.013 


4 


10.5 ± 1.6 


16.2 ± 2.5 


50.4 ± 4.6 




28.1 ± 5.5 


28.7 ± 3.8 


49.6 ± 3 


3 


25.9 ± 4.3 


51.2 ± 8.8 


30 


8 ± 5.6 


14.3 ± 3.2 


16.2 ± 2.3 


10.5 ± 2.9 


43.8 


1120.248 


4 


9.9 ± 1.3 


18. B ± 3.0 


36.4 ± 3.3 


11.0 ± 3.6 


39.2 ± 6.1 


29.1 ±4.6 


43.8 ± 3 


3 


21.0 ± 2.7 


59.1 ±5.7 


62 


4 ± 5.5 


17.8 ±4.1 


11.7 ± 2.3 


5.6 ± 3.0 


36.4 


1120.399 


5 






36.4 ± 3.5 








23.9 ± 4 


2 


6.9 ± 2.6 


16.7 ± 4.6 


33 


5 ± 5.6 


3.2 ± 2.6 


3.7 ± 2.0 




12.0 


1125.540 


5 



































Table 4. Column Density Results 



Sightlii 



N(J = 0) 



N(J = 1) 



HD 24534 


20 


76 


± 





03 


20 


42 


± 





06 


HD 27778 


20 


64 


± 





05 


20 


27 


± 





10 


HD 38087 


20 


39 


± 





08 


20 


29 


± 





05 


HD 40893 


20 


27 


± 





05 


20 


28 


± 





05 


HD 46056 


20 


40 


± 





06 


20 


35 


± 





06 


HD 46202 


20 


38 


± 





07 


20 


38 


± 





07 


HD 63367 


20 


89 


± 





04 


20 


52 


± 





07 


HD 147888 


20 


39 


± 





04 


19 


71 


± 





10 


HD 149404 


20 


60 


± 





03 


20 


34 


± 





05 


HD 170740 


20 


60 


± 





05 


20 


52 


± 





11 


HD 179406 


20 


65 


± 





07 


20 


26 


± 





08 


HD 185418 


20 


34 


± 





04 


20 


56 


± 





06 


HD 186994 


19 


18 


± 





06 


19 


37 


± 





03 


HD 192639 


20 


28 


± 





05 


20 


48 


± 





06 


HD 195965 


19 


90 


± 





03 


20 


18 


± 





03 


HD 197512 


20 


27 


± 





05 


20 


44 


± 





06 


HD 199579 


20 


28 


± 





03 


20 


17 


± 





03 


HD 203938 


20 


72 


± 





05 


20 


68 


± 





08 


HD 206267 


20 


64 


± 





03 


20 


45 


± 





05 


HD 207198 


20 


61 


± 





03 


20 


44 


± 





04 


HD 207538 


20 


64 


± 





07 


20 


58 


± 





05 


HD 210839 


20 


57 


± 





04 


20 


50 


± 





04 



JV(,/ = 2) 

18.40 ± 0.01 
18.47 ± 0.02 
18.98 ± 0.02 
18.01 ± 0.05 
18.28tS:°^ 
18.40lS;°3 
19.22 ± 0.02 
18.511S.S? 

18.86 ± 0.02 

18.32lS:°^ 
17.63 ± 0.02 
18 44+° °^ 

IS 1CI+0-03 
18.19_p 

18.08_o „g 

18.98_Q 04 
,„ -,,+0-02 

'-^■'-'--a-03 

IS OK + O-OS 
18.26„ 

18.56 ± 0.02 
,o -,(. + 0.05 



-;v(.7 = 3) 



17.07 ± 0.03 

,7 r,,:+0.04 
1S 94 + 0-03 
17 34+0-06 
17 40+""'' 
17 66+°"^ 
17.66_o 05 
17.11 ± 0.06 
17.13 ± 0.06 

17 73+0 03 
^'■'•'-0.02 

16.68 ± 0.11 
17 26+"-"^ 
17.12"'" 



17.12"f 



,+0.02 
-0.04 
17 52+° °'' 

17 20+° °^ 

17.08 ± 0.07 

9 + 0-01 
-0-02 
17.68 ± 0.07 

16 76+0 0* 

ID. (D_Q 03 

17 00+°"°"* 

i7-62lS:S'3 

16.93tg:°ol 



N(J = 4) 

15 20+° °'^ 

15 64+" "° 

17.25 ± 0.05 

15.21 ± 0.03 
, 7^+0.05 
la. / j_0.06 

15.72 ± 0.03 

i6.25l°o;°o^ 

15.66l°o:°°7 
16.05 ± 0.04 
17.08 ± 0.05 
16.23 ± 0.07 
16.34 ± 0.03 

15.02i°o°^3 

15 7Q+° °3 
10. '»_o.o4 

16.45 ± 0.02 
16.63 ± 0.06 

15.87 ± 0.02 

IK ns+° °** 
16.08_o.o7 

15.39 ± 0.03 

15.73 ± 0.03 

16 04+° °" 
iu.ui_o.06 

3+0.04 
-0.05 



JV(J = 5) 



JV(J = 6) 



N(J = 7) 



± 0.07 
3+0.17 

-0.19 
„+0.10 

-0-09 
, + 0-05 
^-0-04 
f, + 0-04 

-0.07 
^ + 0.04 

-0.05 
tr+O.OS 

-0.09 
„+0.09 

-0.07 
I ± 0.04 



14.21 : 

14. 33:^ 

16.961 

14.521 

15.061 

15.141 

16.351 

16.131 

15.49 : 

15.91 ± 0.17 

14 44+0 07 
14.^4_0 08 
14 sti+0 05 
14.bS_0 04 

14 31+0 05 
14.01_o 03 

15 os+0 05 
ia.^8_o 04 

14.57 ± 0.02 

15 06+° °^ 
10.UD_o 06 
-1 r ^ft + 0-03 
15.46lo 02 
15.52 ± 0.09 

14.91 ± 0.03 

14.89+ 



16.09^ 



Q+0.05 
-0-04 
14.82 ± 0.07 
4-0.04 
-0.03 



15 46+°-^° 14 89+°-^^ 

iD.^D_o.i5 i^-oy_o.ii 



+ 0.39 
-^-1.12 



I4.08l°;i« 

14.24i°o:i| 
14.17 ± 0.04 

14.72l°o:f3 
l3.78±°o;l| 



14 06+°'" 
i«.UD_o 13 

-I, o-,+0.08 

13.81_o.09 

14.38_o.26 
14 31 + 02 
1,1.01_0.03 
14.17l°oi^ 
13.96t°oi? 



+0.14 
-0.18 



N(J = I 



I3.80tg;i° 



15.58"' 



,+0.02 
■^-0.03 



14.18 ± 0.05 



14.29 ± 0.03 



2.6 ± 0.1 
3 n+0.2 

2.4 ± 0.1 
9.6 ± 0.2 

R , + 0.1 

7.2 ± 0.1 
4.6 ± 0.1 

3.1 ± 0.1 

8.3 ± 0.1 

2 4+0-2 
•^■*-o.i 

6.3 ± 0.2 

4.2 ± 0.1 
6.0 ± 0.1 

6 5+°-l 
°°-0.2 

6.4 ± 0.1 

7 2+°-l 
'■''-0.2 

12.5 ± 0.1 

5.3 ± 0.2 
6.6 ± 0.1 
5.6 ± 0.1 

4.6 ± 0.1 

9.7 ± 0.1 



o 
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Table 5. Limiting Column Densities 



Sightline First Undetected J-level J = 0, u = I 





J 


Line (A) 




(mA)^ log 


;Arb 




(mA)='''= 


log 




HD 24534 


6 


1021, 


.21 


2. 


,85 


13, 


.68 


1, 


.15 


12, 


.48 


HD 27778 


6 


1054, 


,00 


2. 


,29 


13, 


.73 


3, 


.18 


12, 


.98 


HD 38087 


9 


1034, 


.29 


3. 


,37 


13, 


.78 




Detected'* 




HD 40893 


6 


1021, 


.21 


3. 


,32 


13, 


.68 


2, 


.35 


12, 


.78 


HD 46056 


7 


1028, 


.80 


7. 


,07 


14, 


.09 


3, 


.60 


12, 


.99 


HD 46202 


6 


1054, 


.00 


3. 


,30 


13, 


.87 


3, 


.96 


13, 


.03 


HD 53367 


7 


1028, 


.80 


7. 


,26 


14, 


.15 


4, 


.00 


13, 


.06 


HD 147888 


7 


1028, 


.80 


4. 


,66 


13, 


,95 


1, 


,63 


12, 


.64 


HD 149404 


8 


1006, 


.67 


4. 


16 


13, 


,62 


1, 


,83 


12, 


.67 


HD 170740 


7 


1028, 


.80 


6. 


,69 


14. 


,32 


5, 


,22 


13, 


.34 


HD 179406 


7 


1028, 


.80 


8. 


,20 


14. 


17 


2, 


,80 


12, 


.87 


HD 185418 


6 


1054, 


.00 


2. 


,30 


13. 


,71 


2, 


14 


12, 


.76 


HD 186994 


6 


1054, 


.00 


2. 


,00 


13. 


,64 


1, 


,56 


12, 


.60 


HD 192639 


7 


1028, 


.80 


8. 


,46 


14. 


18 


4, 


,38 


13, 


.08 


HD 195965 


7 


1028, 


.80 


0. 


,62 


12. 


,94 


0, 


,43 


12, 


.03 


HD 197512 


7 


1028, 


.80 


9. 


,76 


14. 


,25 


3, 


,92 


13, 


.02 


HD 199579 


8 


1006, 


.67 


1. 


,82 


13. 


,24 




Detected'* 




HD 203938 


7 


1060, 


.03 


6. 


17 


14. 


,21 


8, 


,99 


13, 


.50 


HD 206267 


7 


1005, 


.80 


3. 


,82 


13. 


,77 


2, 


,47 


12, 


.81 


HD 207198 


6 


1021, 


.21 


6. 


,47 


14. 


,05 


3, 


,36 


12, 


.96 


HD 207538 


6 


1054, 


.00 


4. 


,38 


14. 


,03 


3, 


14 


12, 


.94 


HD 210839 


8 


1006, 


.67 


1. 


,47 


13. 


14 


1, 


,32 


12, 


.52 



^This is the l-cr limit on Wx 

3-<j limit on the column density, l-cr limits on W\ are multiplied by 
3 (per Eq. [TJ, and then that upper limit is matched to the curve of growth 
fe- value from Table |4l 

'^Based on a line at 1052.73 A 



<*See Tabled 
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Table 6. Vibrationally Excited Column Densities 



Sightline J v Line (mA) W\ logiV^ 



HD 38087 


1 


1052.73 


16.26 ±5.02 


J-O.0D_Q 25 


1 


1 


1052.61 


13.55 ±5.58 


13.56+[J:3^ 


1 


1 


1097.86 


4.37 ±2.63 


1 q 97 + 0:24 

J->J-^'_0.43 


1 


1 


1128.22 


18.06 ± 5.44 


i4.oot[i:2^ 


Alternate fit'" 








13.54t°i^ 


Alternate fif^ 








13.45;ojJ 


HD 199579 


1 


1052.73 


5.04 ±1.84 


12.62t°i^ 



^ J = 0, = 1 column densities are calculated by fitting the 
equivalent width to the corresponding curve-of-growth with 
the 5- value from the solution in Tabled save for the exceptions 
noted. 

^An alternate solution of the J = 1, = 1 column density 
for HD 38087 based on a curve of growth fit assuming the 
5-value from Table H] for all lines simultaneously instead of 
individually. 

'^An alternate solution of the J = 1, 1/ = 1 column density 
for HD 38087 based on an independent curve of growth fit for 
the three measured lines. The fit for the 6- value is 30 km s~^, 
which should be interpreted to mean that the best fit places 
all lines on the linear portion of the curve of growth (see dis- 
cussion in §3.1.ip . 
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Table 7. Correlations with Other Species — Column Densities^ 



Element 


Ar(H2), total 


Af(2)/Af(0) 


N(3)/N(l) 


Af(4)/Af(0) 


JV(6)/JV(1) 


Af(4)/Af(2) 


JV(5)/JV(3) 


JV(6)/JV(2) 


Li I 


-0.477 


-0.065 


-0.046 


-0.002 


0.217 


0.049 


0.321 


-0.126 




83.7% 


14.2% 


10.1% 


0.6% 


45.3% 


10.7% 


63.4% 


21.1% 




10 


10 


10 


10 


10 


10 


10 


7 


K I 


0.504 


-0.542 


-0.651 


-0.559 


-0.525 


-0.308 


-0.174 


-0.602 




95.4% 


97.0% 


99.4% 


97.6% 


96.3% 


75.4% 


48.2% 


86.1% 




16 


16 


16 


16 


16 


16 


16 


10 


Ca I 


-0.133 


-0.227 


-0.192 


-0.135 


-0.041 


0.016 


0.100 


0.187 




37.7% 


60.3% 


52.5% 


38.2% 


12.0% 


4.8% 


28.7% 


41.7% 




16 


16 


16 


16 


16 


16 


16 


11 


Fe I 


-0.197 


-0.070 


0.072 


-0.400 


-0.022 


-0.470 


-0.087 


-0.431 




43.9% 


16.2% 


16.7% 


77.8% 


5.2% 


85.6% 


20.2% 


71.3% 




11 


11 


11 


11 


11 


11 


11 


8 


Fe II 


-0.108 


0.001 


0.184 


-0.139 


-0.135 


-0.200 


-0.309 


-0.274 




36.8% 


0.4% 


58.7% 


46.3% 


45.0% 


62.7% 


83.8% 


66.7% 




22 


22 


22 


22 


22 


22 


22 


14 


HD 


0.730 




-0.294 


-0.448 


-0.696 


-0.247 


-0.441 


-0.680 




98.9% 


04. ()'/;. 


(il.9% 


83.3% 


94.7% 


53.6% 


82.5% 


77.2% 




11 


11 


11 


11 


11 


11 


11 


6 


CH 


0.870 


-0.291 


-0.456 


-0.354 


-0.271 


-0.261 


0.039 


-0.801 




100.0% 


77.3% 


95.0% 


86.3% 


73.8% 


71.9% 


12.5% 


99.8% 




19 


19 


19 


19 


19 


19 


19 


12 


CH+ 


0.285 


0.266 


0.108 


0.471 


0.489 


0.448 


0.518 


0.097 




76.4% 


73.0% 


34.1% 


95.8% 


96.6% 


94.5% 


97.7% 


23.6% 




19 


19 


19 


19 


19 


19 


19 


12 




0.195 


-0.302 


0.056 


-0.760 


-0.558 


-0.641 


-0.661 


-0.603 




41.1% 


60.4% 


12.3% 


98.9% 


90.6% 


95.4% 


90.1% 


61.2% 




10 


10 


10 


10 


10 


10 


10 


6 


C2 


0.525 


-0.278 


-0.431 


-0.353 


-0.621 


-0.232 


-0.211 


-0.294 




90.3% 


59.3% 


81.4% 


71.2% 


89.9% 


50.8% 


46.6% 


47.8% 




11 


11 


11 


11 


11 


11 


11 


7 


CN 


0.455 


-0.262 


-0.065 


-0.137 


-0.161 


0.039 


-0.121 


-0.071 




93.4% 


69.1% 


19.6% 


39.9% 


46.3% 


11.9% 


35.6% 


16.4% 




17 


17 


17 


17 


17 


17 


17 


11 



^For each element there are three rows. The first row is the Pearson correlation coefficient r. The second row is the confidence 
level that r ^ for a f-distribution of the quantity r(n — 2)^'^^(1 — r^)"^*^^ for n — 2 degrees of freedom, where n is the number of 
data points. The third row is n. Combinations with fewer than three available data points are omitted. 
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Table 8. Correlations with Other Species — Other Elemental Ratios' 



Element 


]V(H2), total 


iV{2)/JV(0) 


JV(3)/JV(1) 


JV(4)/JV(0) 


JV(5)/JV(1) 


]V(4)/]V(2) 


N(5)/N(3) 


]V(6)/]V(2) 


CH+/CH 


-0.311 


0.433 


0.397 


0.660 


0.621 


0.577 


0.445 


0.461 




80.5% 


93.6% 


90.8% 


99.8% 


99.5% 


99.0% 


94.4% 


86.9% 




19 


19 


19 


19 


19 


19 


19 


12 


CH+/CN 


-0.321 


0.346 


0.132 


0.381 


0.365 


0.193 


0.284 


0.002 




79.1% 


82.6% 


38.7% 


86.8% 


85.0% 


54.2% 


73.1% 


0.5% 




17 


17 


17 


17 


17 


17 


17 


11 


Fc II/Fc I 


0.148 


0.167 


-0.074 


0.401 


0.109 


0.397 


0.203 


0.353 




33.6% 


37.5% 


17.1% 


77.8% 


25.1% 


77.3% 


45.0% 


61.0% 




11 


11 


11 


11 


11 


11 


11 


8 



''^Sec notes on Table [7] 
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Table 9. Correlations with Other Species — J Levels' 



Element 


J = 


J = 1 


J = 2 


J = 3 


J = 4 


J = 5 


J = 6 


J = 7 


Li I 


-0.531 


-0.033 


-0.612 


-0.369 


-0.206 


-0.120 


-0.768 






85.9% 


6.7% 


92.0% 


67.2% 


40.5% 


24.1% 


92.5% 






10 


10 


10 


10 


10 


10 


7 




K I 


0.430 


0.474 


-0.444 


-0.550 


-0.384 


-0.441 


-0.732 


-0.929 




89.1% 


92.6% 


90.3% 


96.7% 


84.2% 


90.0% 


97.5% 


75.9% 




16 


16 


16 


16 


16 


16 


10 


4 


Ca I 


-0.377 


0.157 


-0.588 


-0.268 


-0.292 


-0.068 


-0.347 






83.4% 


42.3% 


97.9% 


66.7% 


70.9% 


18.9% 


67.3% 






16 


16 


16 


16 


16 


16 


11 




Fe I 


-0.062 


-0.621 


-0.707 


-0.760 


-0.715 


-0.747 


-0.604 






13.5% 


94.4% 


97.8% 


98.9% 


98.0% 


98.7% 


84.9% 






11 


11 


11 


11 


11 


11 


8 




Fe II 


-0.344 


-0.250 


-0.201 


0.105 


-0.224 


-0.217 


-0.120 


0.260 




87.3% 


72.6% 


61.7% 


34.9% 


67.0% 


65.6% 


30.3% 


16.8% 




22 


22 


22 


22 


22 


22 


14 


4 


HD 


0.691 


0.566 


0.119 


0.056 


-0.143 


-0.591 


-0.811 






97.3% 


91.2% 


25.6% 


12.3% 


30.7% 


92.8% 


90.4% 






11 


11 


11 


11 


11 


11 


6 




CH 


0.796 


0.614 


0.354 


0.057 


0.178 


0.018 


-0.326 


-0.984 




100.0% 


99.3% 


85.1% 


17.7% 


52.1% 


5.8% 


67.2% 


88.7% 




19 


19 


19 


19 


19 


19 


12 


4 


CH+ 


-0.045 


0.321 


0.281 


0.467 


0.583 


0.700 


0.391 


0.564 




14.0% 


80.6% 


74.2% 


94.9% 


98.9% 


99.9% 


76.5% 


38.2% 




19 


19 


19 


19 


19 


19 


12 


4 


l^CO 


0.601 


-0.312 


0.058 


-0.107 


-0.624 


-0.887 


-0.896 






91.3% 


58.7% 


11.8% 


21.7% 


92.7% 


99.9% 


89.6% 






10 


10 


10 


10 


10 


10 


5 




C2 


0.437 


0.623 


-0.174 


-0.081 


-0.273 


-0.322 


-0.737 






79.3% 


94.6% 


37.0% 


17.5% 


55.5% 


63.6% 


90.5% 






11 


11 


11 


11 


11 


11 


7 




CN 


0.769 


0.006 


0.197 


-0.050 


0.159 


-0.160 


0.192 






100.0% 


1.8% 


53.6% 


14.6% 


44.4% 


44.7% 


40.5% 






17 


17 


17 


17 


17 


17 


11 





^Scc notes on Tabic [71 The correlation coefficients in this tabic are partial correlation 
coefficients that assume iV(Htot) is held fixed to reduce the effect of spurious correlations 
arising from total column density effects. 



- 56 - 



Table 10. Excitation Temperatures 



Sightline Toi^ To2^ 



HD 24534 
HD 27778 
HD 38087 
HD 40893 
HD 46056 
HD 46202 
HD 53367 
HD 147888 
HD 149404 
HD 170740 
HD 179406 
HD 185418 
HD 186994 
HD 192639 
HD 195965 
HD 197512 
HD 199579 
HD 203938 
HD 206267 
HD 207198 
HD 207538 
HD 210839 



57 ± 3 73 ± 

55 ± 5 78 ± 

70 ± 6 120 ± 
78 ± 6 74 ± 

73 ± 6 79 ± 
77 ± 8 83 ± 
55 ± 3 96 ± 
45 ± 3 88 ± 
60 ± 3 72 ± 

71 ±8 91 ± 
59 ± 5 67 ± 
100 ±9 79 ± 

96 ± 8 93 ± 

97 ±9 85 ± 
109 ±7 90 ± 
94 ± 8 78 ± 
69 ± 3 76 ± 

74 ± 7 91 ± 

64 ± 3 68 ± 

65 ± 3 73 ± 

72 ± 6 81 ± 
72 ± 4 70 ± 



1 


211 ± 5 


1 


192 ± 7 


4 


460 ± 11 


2 


237 ± 6 


2 


263 ± 8 


2 


247 ± 6 


2 


282 ± 10 


1 


340 ± 13 


1 


406 ±9 


2 


509 ± 34 


2 


342 ± 17 


1 


220 ±4 


2 


215 ±3 


2 


300 ±7 


2 


262 ±3 


2 


328 ± 13 


1 


439 ±4 


2 


302 ± 13 


1 


336 ±7 


1 


284 ±8 


2 


225 ±5 


1 


545 ±8 



'^Temperature derived through hnear fit of 
Io^ NjIqj for J = — 1, column densities 
from iRachford et al] (|2002bl . bood ) . 

^Temperature derived through linear fit of 
logA^j/5J for J = 0-2. 



'^Temperature derived through linear fit of 
\ogNj/gj for J > 3. 
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Table 11. Model Input Parameters' 



Parameter HD 46056 HD 147888^ 



Av 


1.30 


1.94 


Eb-v 


0.49 


0.47 


Rv 


2.66 


4.13 


yo 


4.581 


4.587 


7 


0.911 


1.022 


Cl 


-0.641 


1.611 


C2 


0.878 


0.133 


C3 


3.559 


3.823 


C4 


0.609 


0.339 


iV(Htot) 


3.4 X 10^1 


5.7 X 10 


^tot)/ Eb-v 


6.9 X 10^1 


1.2 X 10 



^Extinction c urve parameters are from 
fl2004l l A^CHtnt.) is from 



Valencic et al 



Lyma n-o; measurements: iDiplas &: Savage 



(1994) for HD 46056 and ICartledge et al. 
( 120041 ) for HD 147888 (note the extensive 
discussion in that paper due to concerns 
about stellar contamination). 

"^Slight discrepancies exist for Ay 
and Rv for HD 147888 between this 
table and Table [2] due to differences 



in 


Rachford et al. 


( 


2002fc 




2009 


) and 


Valencic et al 




2004 


) 


Rachford et al. 


f 2002b 




2009 


) values were preferred for 



Table El but the IValencic et al.l (120041 ) 



values were the published values available 
at the time the models were calculated. 
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Table 12. Modeling Output Parameters 





HD 46056 


HD 


147888 


Parameter 


Observations^ 


Best Model Fit 


Observations'* 


Best Model Pit 


N(0) 


20.40 ± 0.06 


20.99 


20.39 ± 0.04 


21.28 


N(l) 


20.35 ± 0.06 


20.80 


19.71 ± 0.10 


20.91 


N(2) 


l8.28t°;S5 


18.88 




18.45 


N(3) 




17.65 


17.11 ± 0.05 


16.68 


iV(4) 




15.89 


15 65+° '"^ 


15.72 


iV(5) 


15 06+° '^* 


15.48 


15 13 + 009 


15.26 


iV(6) 


14 11 + 0.39 


14.52 


14 24+° '-'' 


14.33 


log JV(CH) 


13.31 ± 0.06 


14.28 


13.30 ± 0.02 


14.59 


log JV(CH + ) 


13.37 ± 0.05 


11.50 


12.91 ± 0.04 


11.92 


log JV(12C0) 




14.48 


15.30 ± 0.06 


14.43 


log Ar(CN) 


12.30 ± 0.10 


12.36 


12.35 ± 0.02 


12.20 


/{H2) 


0.285 ± 0.087 


0.959 


0.105 ± 0.025 


0.974 


JV(H2°)/iV(H2P) 


0.885 ± 0.012 


0.640 


0.209 ± 0.006 


0.430 


nn ( cm~^) 




5 X 10^ 




3 X 10^ 


Radiation (G) 




3 X 10^ 




1 X 10^ 


TcORB (K) 




46 




44 



^Observed column densities and subsequently derived quantities. H2 from this paper, other 
column densities from http: //www .astro. uiuc.edu/ ^dwelty/coldens.htmlj 



